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Abstract 


Leukocyte  (WBC)  adhesion  to  venular  endothelium  (EC)  is  an  important 
step  in  the  immune  response  which  precedes  their  emigration  through  the 
microvessel  wall.  To  assess  the  relationship  between  membrane  WBC 
mechanics  and  the  formation  and  disruption  of  WBC-EC  bonds, 
measurements  of  transients  in  contact  area  and  WBC  shape  were  made 
during  adhesion  in  postcapillary  venules  of  the  rat  mesentery  during  control 
conditions  (spontaneous  adhesion),  in  response  to  exposure  to  the 
chemoattractant  N-formyl-methionyl-leucyl-phenylalanine  (FMLP)  and 
following  tissue  exposure  to  interleukin- 1/9  (iL-1).  To  assess  the  roles  of 
WBC  deformability,  adhesion  was  also  studied  during  suffusion  with 
colchicine  and  cytochalasin  B  (CB).  Frame-by-frame  analysis  of  video 
recordings  facilitated  measurement  of  the  length  of  the  contact  zone 
between  WBC  and  EC  (L^)  and  the  acute  angle  (9)  between  the  trailing  edge 
of  the  WBC  membrane  and  EC  surface.  During  the  initial  adhesion,  was 
found  to  steadily  increase  with  time  as  new  bonds  were  formed,  and  to 
subsequently  decrease  with  time  as  the  trailing  edge  of  the  WBC  membrane 
began  to  peel  away  from  the  EC.  The  contact  angle  at  the  trailing  edge  of 
the  cell  decreased  exponentially  for  control,  IL-1,  CB,  and  colchicine.  For 
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FMLP  the  angle  initially  decreased  then  subsequently  increased  as  the  WBC 
became  more  activated.  Both  6  and  plotted  against  furnished  a 
measure  of  deformabiiity  for  control,  FMLP,  IL-1,  and  colchicine  which  was 
consistent  with  current  data.  6  decreased  with  for  control,  IL-1,  and 
colchicine  treatment  indicating  the  cells  were  deformable.  0  increased  with 
for  FMLP  indicating  that  the  more  activated  and  stiffest  ceils  adhered  at 
the  highest  shear  rates.  L^  versus  provided  consistent  results  with  the 
exception  of  IL-1  stimulated  adhesion  in  which  L^  was  invariant  with  r^.  Tiie 
time  of  transition  between  bonding  and  peeling,  tp,  was  assumed  to 
correspond  to  the  iife  span  of  bonds  formed  during  the  initial  WBC-EC 
contact.  A  logarithmic  decline  of  tp  with  wall  shear  stress,  r^,  suggested  a 
mode  of  peeling  consistent  with  the  kinetic  theory  of  fracture.  The 
logarithmic  relationship  was  significant  for  peeling  under  control,  FMLP,  and 
IL-1  conditions.  Employing  the  kinetic  theory  of  fracture,  the  bond  energy 
(U„)  and  bond  force  (f^)  were  calculated  for  each  condition  and  were  3% 
greater  for  FMLP  as  compared  to  control  and  IL-1  (p  <  0.05).  U„  for  both 
control  and  IL-1  was  not  significantly  different  (p  >  0.05),  indicating  that 
both  stimuiants  for  adhesion  induce  WBC  sticking  through  similar 
mechanisms.  To  model  the  adhesion  process  and  elucidate  the  relative  roles 
of  cell  deformabiiity  and  strength  of  the  adhesion  bond  as  determinants  of 
the  peeiing  process,  a  theoretical  analysis  was  performed  in  which  the  bond, 
normal  stress,  and  membrane  tension  distributions  were  anaiyzed.  Since  the 
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avidity  of  the  receptor-ligand  bonds  for  FMLP  and  IL-1  adhesion  were  similar 
in  magnitude,  the  FMLP  adhesion  appeared  to  be  much  stronger  due  to 
increased  bond  density  rather  than  bond  force. 
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Chapter  1 
Introduction 


1.1  Rationale 

Adhesion  of  leukocytes  to  post-capillary  venular  endothelium  Is  a 
central  component  of  the  immune  response  which  precedes  the 
extravasation  of  white  blood  cells  (WBCs)  through  the  vessel  wall. 
Elucidation  of  the  mechanics  of  the  adhesion  process  is  necessary  to 
understand  the  etiology  of  pathological  conditions  which  result  from  either 
weak  or  excessive  WBC  adhesion  to  the  endothelium. 

For  example,  leukocyte  adhesion  deficiency  syndrome  results  from  the 
absence  of  a  specific  adhesion  molecule  and  manifests  itself  in  virulent 
bacterial  infections,  and  early  death  (Lasky,  1992).  Excessive  leukocyte 
adhesion  may  also  have  deleterious  consequences.  Progressive  leukocyte 
adhesion  has  been  implicated  In  the  "no  re-flow  phenomenon"  following 
periods  of  myocardial,  brain,  and  bowel  ischemia  and  acute  hemorrhagic 
shock  (Bagge,  1984;  Barroso-Arnanda  ef  a/.,  1988;  Dahlgren  ef  a/.,  1984; 
Engler  et  a/.,  1983  and  1986;  Engler,  1989,  Granger  et  aL,  1989). 
Concomitant  to  the  period  of  ischemia,  there  is  a  progressive  increase  in 
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vascular  resistance  due  to  the  continual  entrapment  of  WBCs  in  the 
microcirculation  and  edema  formation  (Engler  et  a!.,  1 986).  The  entrapment 
or  leukocyte  capillary  plugging  occurs  at  points  within  the  capillary  bed  that 
require  the  WBC  to  deform  in  order  to  pass,  such  as,  capillary  entrances  and 
tapered  segments  caused  by  protruding  endothelial  cell  nuclei  (Dahlgren  et 
a/.,  1 984).  Upon  reperfusion,  hemodynamic  forces  are  insufficient  to  disrupt 
the  adhesion  between  the  leukocyte  and  endothelium;  therefore,  a  "no- 
reflow  phenomenon"  occurs  and  tissue  ischemia  is  exacerbated.  While 
trapped  in  the  microcirculation  and  parenchymal  tissue,  the  activated  WBCs 
may  release  a  number  of  toxic  substances,  such  as,  lytic  enzymes  and 
oxygen  free  radicals,  resulting  in  endothelical  injury  and  edema  which  further 
increases  vascular  resistance  (Barroso-Arnanda  et  a!.,  1 988).  Anti-adhesive 
treatments  are  currently  being  investigated  to  aid  in  the  prevention  of  the 
"no-reflow  phenomenon"  and  reduction  of  tissue  injury  following  ischemic 
episodes  (Engler  et  a!.,  1986;  Suzuki  et  a!.,  1991;  Granger  et  at.,  1989; 
Schmid-Schonbein  and  Engler,  1990).  Leukocyte  mediated  injury  has  also 
been  implicated  in  allograft  rejection,  rheumatoid  arthritis,  inflammatory  skin 
disease,  and  adult  respiratory  distress  syndrome  (Harlan  et  a/.,  1992).  For 
these  reasons,  the  study  of  WBC  adhesion  is  a  pertinent  and  compelling 
topic  of  microvascular  research. 

In  vivo,  leukocytes  can  be  observed  rolling  along  the  vascular 
endothelium,  particularly  the  venules.  During  inflammation,  the  WBCs 
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adhere  to  the  endothelium  and  emigrate  across  the  microvessel  wall  to  the 
site  of  injury  (Grant,  1973).  Their  adhesion,  and  ultimately  the  progression 
of  the  inflammatory  response,  rests  upon  the  balance  of  the  adhesive  and 
hemodynamic  forces  acting  upon  the  leukocyte.  Various  agents  may  be 
employed  to  alter  adhesive  and  hemodynamic  forces  acting  on  the  WBC. 
There  are  several  modes  of  adhesion,  two  of  which  can  be  stimulated  with 
separate  agents,  N-formyl-methionyi-ieucyi-phenylalanine  (FMLP)  and 
interleukin-1/9  (IL-1 );  therefore,  possibly  manifesting  distinct  adhesive  forces. 
Hemodynamic  forces  can  be  altered  by  affecting  WBC  deformability.  Two 
agents  are  specifically  proposed  for  this  purpose;  cytochalasin-B  and 
colchicine,  both  which  cause  an  increase  in  cell  deformability  by  acting  on 
the  WBC's  cytoskeleton.  In  addition  FMLP  has  been  shown  to  significantly 
decrease  WBC  deformability.  In  studying  the  relationship  between  adhesive 
and  hemodynamic  forces  under  these  different  conditions,  a  greater 
understanding  of  the  mechanical  properties  of  the  adhesion  molecules  may 
provide  new  insight  into  the  therapeutic  management  of  leukocyte- 
endothelium  adhesion. 

1.2  Background 

Historically,  the  strength  of  the  adhesive  force  between  leukocytes 
and  vascular  endothelium  has  been  estimated  indirectly  by  measurements  of 
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white  blood  cell  rolling  velocity  (V^)  and/or  the  number  of  cells  per  unit  area 
of  endothelium.  Atherton  and  Born  (1973)  concluded  that  leukocyte  rolling 
was  governed  by  two  forces:  the  hemodynamic  shear  force  and  the  adhesive 
force  between  WBC  and  endothelium.  Further,  if  the  adhesive  force  was 
assumed  to  be  similar  for  all  WBCs,  the  proportionality  between  the  velocity 
of  blood  (red  blood  cell  velocity,  V^h,)  and  was  due  to  the  balance  of 
these  shear  forces.  Mayrovitz  and  Wiedeman  (1976)  showed  that  wall  shear 
stress  was  indeed  an  important  parameter  in  determining  whether  leukocyte 
adherence  would  occur.  Using  leukocyte  flux  in  arterioles  (bat  wing)  as  a 
measure  of  the  strength  of  adhesion,  they  concluded  that  the  initial 
adherence  of  WBCs  was  dependent  upon  the  magnitude  of  shear  stress,  up 
to  a  critical  or  threshold  value  of  8  dyne/cm*,  above  which  no  WBC 
adhesion  could  be  observed.  Contrastingly,  they  found  that  once  the  cells 
were  rolling  along  the  endothelium,  there  was  no  clear  relationship  between 
V„bc  3nd  shear  stress  except  for  at  low  blood  velocities.  Similar  findings 
were  published  by  Firrell  and  Lipowsky  (1989)  based  upon  observations  in 
postcapillary  venules  (mesentery).  The  invarience  of  V^bc  with  shear  stress 
indicated  that  further  mechanisms  needed  to  be  considered,  such  as,  cell  to 
cell  adhesive  forces  and  cell  deformability. 

Cell  to  cell  adhesive  forces  are  derived  from  the  interaction  of 
receptors  on  the  WBC  and  endothelial  cells.  Preliminary  evidence  indicates 
that  "rolling  receptors"  or  selectins  enable  neutrophils  to  roll  along  the 
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endothelium.  The  slowing  of  the  white  cell  and  the  close  proximity  to  the 
vascular  wall,  in  turn,  promotes  the  engagement  of  the  integrins  (adhesion 
molecules),  thereby,  facilitating  adhesion  and  ultimately  cell  emigration, 
Lawrence  and  Springer  (1 991)  and  Figdor  and  Van  Kooyk  (1992).  Leukocyte 
deformability  may  aid  in  the  adhesion  process  by  helping  to  balance  the 
adhesive  and  hemodynamic  forces.  The  more  deformable  the  WBC,  the 
longer  the  cell  will  stretch  and  the  larger  the  area  of  contact  between  the  cell 
membrane  and  endothelium.  Contact  area  may  be  the  principal  determinant 
of  the  maximal  force  that  can  be  withstood  before  hemodynamic  stresses 
disrupt  the  bond  between  WBC  and  EC.  Schmid-Schonbein  et  aL  (1987) 
presented  an  analysis  affirming  this  relationship.  They  showed  that  the  force 
per  unit  length  it  takes  to  break  the  adhesive  bonds  (fracture  stress)  between 
a  steady  rolling  white  cell  and  endothelium  is  inversely  proportional  to  the 
ratio  of  Therefore,  the  lower  the  velocity  ratio;  the  higher  will  be 

the  fracture  stress.  Since,  the  fracture  stress  (adhesion  energy  density)  is 
equal  to  adhesion  energy  divided  by  contact  area,  a  larger  contact  area,  will 
yield  a  greater  adhesive  energy  between  the  WBC  and  EC  for  a  given 
fracture  stress.  Therefore,  a  highly  deformable  cell  may  have  a  larger 
contact  area  which  will  increase  the  strength  of  the  cell  to  cell  adhesion. 
Firrell  and  Lipowsky  (1989)  hypothesized  that  the  stretching  and  elongation 
of  the  WBC  in  shear  flow  may  also  help  temper  hemodynamic  shear  forces 
by  reducing  the  extent  of  the  vessel  lumen  obstructed  by  the  cell;  thereby. 
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reducing  the  wall  shear  stress  acting  on  the  ceil.  Thus,  cell  deformability 
influences  WBC-endothelium  interaction  by  increasing  contact  length  and 
decreasing  hemodynamic  shear  forces,  hence,  influencing  the  balance  of 
bond  and  hemodynamic  forces  in  an  adherent  ceil. 

It  can  be  seen  by  the  research  discussed  above  that  contact  area  is  an 
important  determinant  in  assessing  adhesive  strength,  yet  little  attention  has 
been  given  to  this  parameter.  Schmid-Schonbein  et  a/.  (1975)  estimated 
contact  lengths  of  adhered  ceils  in  venules  of  the  rabbit  omentum  in  order 
to  estimate  the  shear  force  on  the  WBCs  due  to  differing  blood  hematocrits. 
Firrell  and  Lipowsky  made  measurements  of  contact  length  {LJi  as  a  function 
of  wall  shear  rate  finding  that  increased  in  an  exponential  fashion  as  wall 
shear  rate  increased.  Yet,  there  are  no  published  reports  to  date  examining 
contact  length  change  over  time  in  an  adhered  cell. 

The  aim  of  this  present  experimental  study  was  to  look  specifically  at 
the  transients  in  contact  length  (LJ  between  the  WBC  and  endothelium  and 
acute  angle  of  contact  {0)  the  WBC  ♦orms  with  the  endothelium  during 
spontaneous  intermittent  adhesion  and  relate  these  parameters  to  bond  life 
time  and  cell  deformability.  It  was  hypothesized  that  the  breaking  of  the 
bonds  responsible  for  adhesion  in  shear  flow  would  follow  the  kinetic  theory 
of  fracture.  This  relationship  served  to  elucidate  the  mechanical  strength 
parameters  of  the  adhesion  bonds  which  in  turn  were  used  as  part  of  a 
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theoretical  analysis  of  WBC  membrane  mechanics  assuming  force  equilibrium 
between  hemodynamic  and  adhesive  forces. 


Chapter  2 

Experimental  Methods  and  Data  Acquisition 
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2. 1  Methods 

Female  Wistar  rats,  100-200  g  in  weight  were  anesthetized 
intraperitoneally  with  a  Nembutal  pentobarbital  sodium  (Abbott  Laboratories) 
injection  (35  mg/kg).  After  a  tracheostomy,  the  left  jugular  vein  and  right 
carotid  artery  were  cannulated  with  polyethylene  tubing  (PE-50)  for 
administration  of  supplemental  anesthetic  (1/1 0th  of  original  dose)  and 
monitoring  of  blood  pressure,  respectively. 

A  Century  Technology  CP-02  Physiologic  Pressure  Transducer  was 
connected  to  the  rats  carotic  artery  catheter  and  to  a  Beckman  Instruments 
Dynograph  Recorder  R611.  Recorder  components  included  TYPE  9853A 
Voltage/Pulse/Pressure  Coupler,  preamplifier  type  461 D,  and  power  amplifier 
type  R-41 1 .  Blood  pressure  was  monitored  throughout  the  experiment  to 
assess  physiologic  state  and  level  of  anesthesia. 

A  midline  abdominal  incision  was  made  and  a  section  of  the  distal 
intestinal  mesentery  exteriorized.  The  rat  was  placed  on  a  plexiglass  stage, 
and  the  intestine  and  associated  mesentery  were  draped  on  a  glass  slide  for 
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direct  observation  with  transmitted  light.  Figure  2.1  is  a  schematic  of  the 
experimental  set-up.  The  intestine  was  handled  gingerly  to  prevent  not  only 
tissue  damage  but  peristalsis  due  to  over  stimulation.  The  tissue  was 
continuously  suffused  with  warmed  Ringer's  gelatin  solution  (37°  C)  buffered 
with  sodium  bicarbonate  and  adjusted  to  a  pH  of  7.4  with  hydrochloric  acid. 
The  mesentery  was  allowed  to  stabilize  for  a  20  minute  period  prior  to  the 
start  of  data  acquisition.  If  an  abundance  of  leukocyte  sticking  or  significant 
stasis  in  the  microvessels  was  observed,  the  mesenteric  section  was  not 
used. 

Venules  were  observed  at  high  magnification,  Zeiss  40/0.75  water- 
immersion  objective  (magnification/numerical  apperature)  and  projected  by 
a  79  mm  eyepiece  onto  a  silicon  diode  video  camera  (DAGE-MTI  SIT-66X). 
The  images  were  then  recorded  on  1/2  inch  VHS  videotape  (Panasonic 
Omnivision  II  VHS  Video  Cassette  Recorder  NV-8950)  for  off-line  analysis. 
A  video  time  code  generator  (FOR. A  Video  Timer  VTG-33)  was  used  to 
provide  a  time  reference  for  off-line  analysis.  White  cell  saltation  and 
intermittent  adhesion  were  recorded  while  focusing  upon  the  lateral  edges 
of  post-capillary  venules.  During  each  adhesion  event,  the  transient 
deformations  of  the  WBCs  were  recorded  for  subsequent  frame  by  fame 
analysis. 

These  observations  were  repeated  for  five  different  protocols:  (1) 
Spontaneous  adhesion  during  a  "control”  state,  (2)  adhesion  in  response  to 


Figure 
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topical  application  of  the  chemoattractant  N-formyl-methionyl-leucyl- 
phenylalanine  (FMLP)  at  a  concentration  of  10'^M  dissolved  in  the  Ringer's 
irrigation  solution,  (3)  adhesion  in  response  to  stimulation  by  the  cytokine, 
interleukin-1/9  (IL-1),  administered  to  the  animal  Intraperitoneally  at  50  U  in 
1  ml  normal  saline  per  hour  for  two  hours  prior  to  mesentery  exteriorization, 
(4)  adhesion  during  suffusion  of  the  tissue  with  Ringer's  solution  containing 
cytochalasin  B  (CB)  at  a  concentration  of  10  ^iM,  (5)  with  the  addition  of 
colchicine  (100//g/mi)  to  the  irrigation  solution. 

Studies  of  adhesion  in  response  to  FMLP  and  IL-1  were  performed  in 
light  of  the  two  different  modes  of  adhesion  stimulated  by  those  agents. 
FMLP  is  a  chemotactic  peptide  that  directly  activates  the  WBC,  and  has  been 
shown  to  significantly  decrease  WBC  deformability  at  concentrations  of  10' 
^M  (Kawaoka  etal.,  1981;  House  and  Lipowsky,  1991;  and  Lipowsky  er  a/. , 
1991).  IL-1  is  a  major  cytokine  that  induces  expression  of  adhesion 
molecules  on  the  endothelium  with  little  effect  on  WBC  adhesion  molecules. 
Resting  endothelium  typically  lacks  the  molecules  to  support  adhesion; 
however,  upon  addition  of  a  cytokine  there  is  an  increase  in  the  expression 
of  three  known  adhesion  molecules;  Intercellular  adhesion  molecule-1  (ICAM- 
1),  vascular  cell  adhesion  molecule-1  (VCAM-1),  and  endothelial  leukocyte 
adhesion  molecule  (ELAM-1)  (Vedas  et  a!.,  1992).  The  deformation 
properties  of  the  WBC  in  response  to  adhesion  due  to  IL-1  endothelial 
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stimulation  have  been  shown  to  decrease;  however,  not  statistically  different 
from  that  of  control  (Lipowsky  et  aL,  1991). 

Cytochalasin  (CB)  and  colchicine  have  both  been  shown  to  increase 
WBC  deformability  (Lipowsky  ef  a/. ,  1 99 1 ) .  Cy tocholasin  B  is  one  of  a  broad 
family  of  cytochalasins  and  is  a  naturally  occurring  organic  molecule  that 
binds  to  actin  and  inhibits  its  polymerization.  It  permeates  the  cell 
membrane  and  causes  the  inhibition  of  polymerization  and  depolymerization; 
therefore,  the  cell  becomes  more  flaccid  (Cooper,  1987).  Colchicine  is  an 
anti-mitotic  agent  commonly  used  In  the  treatment  of  Inflammatory 
disorders;  such  as,  rheumatoid  arthritis.  It  has  been  proposed  that  colchicine 
attenuates  the  inflammation  process  by  disrupting  the  microtubules,  which 
are  a  major  structural  component  of  the  WBC,  thereby,  decreasing  cell 
stiffness.  Increased  WBC  deformability  in  response  to  colchicine  has  been 
confirmed  in  vivo  (Lipowsky  et  al.,  1991  and  Asako  et  ai.,  1992). 

2.2  Data  Analysis 

After  completion  of  the  experiments,  the  video  tape  recordings  were 
analyzed  off-line.  It  should  be  noted  at  this  point  that  several  video 
recordings  were  used  from  previous  experiments  which  utilized  the  same 
protocol.  The  recordings  were  not  previously  analyzed  for  the  information 
desired  in  this  present  study.  Use  of  these  recordings  was  by  permission  of 


the  principal  investigator  of  the  1988  study.  Dr.  H.  H.  Lipowsky.  The 
combination  of  current  and  past  experimental  data  served  to 
provide  a  larger  pool  of  data  and  reduce  the  number  of  animals  required  to 
complete  this  study. 

Intermittently  adherent  cells  in  the  venules,  with  durations  of  adhesion 
which  lasted  for  at  least  2.5  seconds,  but  not  permanently,  were  identified. 
Venule  diameters  was  assessed  by  the  video  image  shearing  technique, 
using  an  IPM  (Instrumentation  for  Physiology  and  Medicine,  Inc)  Image 
Shearing  Monitor  model  908.  Measurements  of  centerline  red  cell  velocity 
(Vrfie)  were  made  using  video  two-slit  photometric  technique.  Using  a 
photoanalyzer  (IPM  Video  Photoanalyzer  model  204),  two  photodetectors 
were  aligned  along  the  venule  of  interest  at  a  spacing  ranging  from  30  to  60 
A/ms.  The  output  of  the  photoanalyzer  was  coupled  to  an  on-line  self¬ 
tracking  correlator  (IPM  Velocity  Tracker  Model  102b).  The  signals  from  the 
photodetectors  were  cross-correlated  to  determine  the  time  shift.  The 
spacing  between  the  detectors  was  then  divided  by  the  time  shift  to  yield 
the 

Estimates  of  mean  blood  velocity  (V^„)  were  then  determined  using  the 
empirical  relationship  peculiar  to  the  two-slit  method  (Lipowsky  and 


Zweifach,  1978). 
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The  wall  shear  rate  (  y  )  estimated  as  the  Newtonian  value, 

y-BV^ID  <2-2) 

where  D  denotes  vessel  diameter.  Wall  shear  stress  (rj  was  calculated  as 
the  product  of  the  wall  shear  rate  and  blood  viscosity  (an  average  value  of 
blood  viscosity  of  0.025  Poise  was  assumed). 

Images  of  WBC  deformation  and  membrane  peeling  were  analyzed  to 
measure  transient  changes  in  white  blood  cell-endothelium  (WBC-EC)  contact 
length  (LJ  and  the  acute  angle  of  contact  {0)  over  the  entire  time  the  cell 
remained  adherent.  Figure  2.2  depicts  both  the  length  and  angle  of  interest. 
The  images  were  digitized  in  successive  frames  by  an  image  sampling 
program,  IMP. EXE,  originally  written  by  Michael  Kent  in  Turbo-C.  The 
program  allowed  the  user  to  sample  up  to  sixteen  consecutive  frames  at  a 
rate  specified  by  the  user,  identify  a  1 28  by  1 28  pixel  region  of  interest 
(ROD  from  each  frame,  and  create  a  montage  of  the  16  ROIs  in  time 
sequence  order.  For  each  cell,  a  minimum  of  sixteen  frames  were  digitized 
over  the  duration  of  the  cell  adhesion  with  more  frames  being  "grabbed"  as 
required  to  adequately  analyze  the  transient  changes  in  and  0. 

The  next  step  in  the  data  analysis  was  to  measure  the  WBC-EC 
contact  length  with  the  measurement  capability  of  the  sampling  program. 


2.2  Schematic  of  adherent  WBC.  9  is  the  contact  angle  between  the  WBC  and  EC.  represents  the  contact 
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However,  the  digitization  process  created  a  distortion  in  the  projected  image 
where  the  number  of  pixels  per  micron  differed  in  the  horizontal  and  vertical 
directions.  Therefore,  the  program  was  first  modified  to  allow  for  separate 
horizontal  and  vertical  calibrations,  thereby,  increasing  the  accuracy  of  the 
contact  length  measurements.  The  separate  calibration  factors  were 
incorporated  into  the  length  calculation  with  the  foiiowing  equation: 

segment  length  =  yJ^HCAL'HLEN)^  +  [VCAL  VLEN)^  (2.3) 

where  HCAL  is  the  horizontal  calibration  factor  (pm/pixel),  HLEN  (pixels)  is 
the  horizontal  length  measurement,  VCAL  is  the  vertical  calibration  factor 
(/ym/pixei),  and  VLEN  (pixels)  is  the  vertical  length  measurement.  At  the 
magnification  used  in  the  present  experiments,  the  width  of  the  entire  video 
fieid  was  approximately  11 4 //m.  Scale  calibration  factors  averaged  0.2278 
//m/pixel  for  HCAL  and  0.1553  //m/pixel  for  VCAL.  The  program  was  also 
enhanced  to  allow  the  automatic  storage  of  iength,  width,  and  time 
measurements  for  each  cell  in  an  ASCII  file  that  could  be  read  directly  into 
a  statistical  program. 

The  contact  angle,  B,  was  measured  by  making  a  thermal  print  of  the 
image  file  and  judging  the  angle  of  the  WBC  membrane  by  eye.  The  angle 
was  measured  the  old  fashion  way,  with  a  protractor.  The  information  was 
then  typed  into  an  analysis  program. 
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2.3  Statistical  Analysis 


Measurements  of  L,,  and  0  were  input  to  Sigma  Plot  scientific  graphing 
system  for  preliminary  analysis  and  graph  generation.  The  data  was  further 
analyzed  using  the  statistical  software  Minitab  Release  7.1. 


Chapter  3 

Experimental  and  Analytical  Results 
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To  assess  the  relationship  between  membrane  WBC  mechanics  and 
the  formation  and  disruption  of  WBC-EC  bonds,  transients  in  contact  area 
and  WBC  shape  were  made  during  spontaneous  adhesion.  Tables  of  the 
hemodynamic  data  for  each  WBC  are  provided  in  Appendix  A  while,  the 
experimental  analysis  for  each  cell  is  listed  in  Appendix  B. 

3.1  Formation  and  Disruption  of  WBC-EC  Bonds 

The  hemodynamic  environment  experienced  by  each  cell  of  interest 
was  very  important  in  assessing  the  adhesive  force  required  to  keep  the 
WBC  adherent  to  the  EC.  Therefore,  the  wall  shear  rate  and  wall  shear 
stress  were  calculated  as  described  in  Chapter  2  for  each  cell  of  interest. 
Table  3.1  lists  the  mean  values  and  standard  deviations  of  the  experimental 
hemodynamic  conditions  for  control,  FMLP,  IL-1,  CB,  and  colchicine 
treatments.  Venule  diameters  ranged  form  16.2  to  43.3  //m.  The  mean 
was  greatest  for  FMLP  and  least  for  colchicine;  indicating  the  WBC  adhesive 
strength  was  greatest  for  the  FMLP  treatment  and  least  for  colchicine. 


Hemodynamic  parameters 
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Typical  examples  of  WBC-EC  contact  length  (L^)  versus  time  for  each 
treatment  are  shown  in  Figure  3.1  (a),  (b),  and  (c)  for  control,  FMLP,  and  IL- 
1  conditions,  respectively.  Figure  3.2  (a)  and  (b)  depict  representative 
examples  of  versus  time  for  treatments  with  CB  and  colchicine, 
respectively.  Since  the  fluid  shear  stress  was  not  controlled  during  the 
experiments,  the  representative  curves  are  at  differing  wall  shear  stresses. 
During  the  initial  adhesion,  steadily  increased  with  time  as  new  bonds 
were  formed  at  the  leading  edge,  and  L,,  subsequently  decreased  with  time 
as  the  trailing  edge  of  the  WBC  membrane  peeled  away  from  the  EC  during 
initiation  of  a  roiling  movement  of  the  WBC.  Linear  regressions  through  both 
the  periods  of  increasing  and  decreasing  enabled  identification  of  the 
maximum  L^.  This  point  was  considered  to  represent  the  equilibrium 
between  bond  formation  and  disruption,  as  the  WBC  began  to  roll  and  the 
corresponding  time  was  denoted  as  the  time  to  peeling,  tp. 

Table  3.2  lists  the  regression  parameters  and  correiation  or  r-value  for 
each  regression  depicted  in  Figures  3.1  and  3.2.  In  general  all  treatments, 
excluding  FMLP,  exhibited  a  greater  rate  of  decrease  in  L,,  than  increase 
which  can  be  interpreted  as  the  bonds  were  broken  faster  than  they  were 
formed.  Figure  3.3  is  a  comparison  of  L^  versus  time  for  control,  FMLP,  and 
IL-1  and  corresponding  tp.  As  demonstrated  in  the  figure,  the  FMLP  treated 
cells  remained  adhered  to  the  endothelium  significantly  longer  than  all  others. 
The  time  to  peeling  is  also  much  greater  for  FMLP  than  control  or  IL-1. 


Figure  3.2  Representative  trends  of  contact  length  versus  time 
for  WBCs  with  (a)  CB  and  (b)  colchicine  treatment. 


vs.  time  regression  parameters 
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and  correspondingly  had  a  greater  t  . 
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Table  3.3  lists  the  mean  and  standard  deviation  of  tp,  at  tp,  and 
total  time  the  cell  remained  adherent  to  the  endothelium  (t,,,)  for  ail  five 
conditions.  As  previously  noted,  r^,  tp  and  t,p,  are  greatest  for  the  FMLP 
treated  cells,  thereby,  indicating  a  strong  bond  adhesive  force. 

Figures  3.4  and  3.5  show  tp  versus  for  all  treatments  on  a  semi-log 
scale.  The  regression  of  In(tp)  and  was  statistically  significant  for  control, 
FMLP,  and  IL-1  for  p  <  0.01.  The  regressions  for  CB  and  colchicine  were 
not  statistically  significant  (p  >  0.05)  due  to  limited  range  of  data. 
Parameters,  t-ratios,  p-values,  and  r-value  are  listed  in  Table  3.4  for  each 
regression. 

Figure  3.6  is  a  composite  of  all  three  significant  regressions;  control, 
FMLP,  and  IL-1 .  The  slopes  for  each  of  the  treatments  were  not  significantly 
different  from  one  another;  while,  the  intercept  for  the  FMLP  treatment  was 
significantly  greater  than  either  the  control  or  IL-1  (p  <  0.05).  The 
logarithmic  relationship  between  tp  and  r*,  suggested  a  mode  of  peeling 
consistent  with  the  kinetic  theory  of  fracture  (Zhurkov,  1965). 

3.1.1  Kinetic  Theory  of  Fracture 

Zhurkov's  theory  considered  the  fracture  of  solids  as  a  time  process 
whose  rate  is  determined  by  mechanical  stress  and  temperature.  A  universal 
rate  relation  was  suggested  as  follows: 


Experimental  time  to  peeling  and  contact  length  measurements 


,  Time  to  peeling  (sec) 


•It 


(a)  Control,  N  =  1 8 


2 

(dyne/cm  ) 

Figure  3.4  Time  to  peeling  versus  well  shear  stress  for  (a)  control 
conditions  and  (b)  FMLP  and  (c)  IL-1  stimulated  adhesion.  All 
regressions  were  statistically  significant  for  p  <  0.01. 


,  Time  to  peeling  (sec) 


Rearession  parameters  and  statistics  for  ln(t„)  versus  r, 
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r  =  r,  exp[ 


{U„-Ya) 

kT 


1 


(3.1) 


where  r  is  the  lifetime  of  a  bond  under  load,  a  is  the  load  tensile  stress,  T 
is  the  absolute  temperature,  k  is  Boltzman's  constant  and  r„  Ug,  and  y  are 
constant  coefficients  which  represent  material  constants  of  the  solid.  For 
the  case  of  constant  temperature,  as  in  the  case  of  the  present  study, 
equation  (3.1)  can  be  rewritten  as: 


ln(,.  =  irnlr.)  .  ^ 


(3.2) 


r,  is  the  reciprocal  of  the  natural  frequency  of  oscillation  of  atoms  in  a  solid. 
Zhurkov  reported  r,  to  be  in  the  order  of  magnitude  of  10'^’  seconds  for  all 
solids,  independent  of  structure  and  chemical  nature.  Using  the  kinetic 
theory  of  strength  of  solids.  Bell  (1978)  calculated  this  parameter  to  be  10'” 
seconds  for  receptor-ligand  bonds.  U,  can  be  interpreted  as  the  magnitude 
of  the  energy  barrier  which  determines  the  probability  of  bond  breakage,  y, 
represents  an  empirical  parameter  that  accounts  for  the  structure  of  the  solid 


and  its  imperfections. 
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The  logarithmic  regression  of  tp  versus  r^, 

In(fp)  =  (3.3) 

can  be  related  to  equation  (3.2).  Therefore,  tp  is  equivalent  to  Zhurkov's  r, 
0^  is  equivalent  to: 

0o  =  ln(r,)  +  ^  (3.4) 

and  0y  is  equivalent  to: 

A  -  ^  •  I3  SI 

C  is  a  constant  relating  the  fluid  shear  stress  to  the  normal  tensile  stress 
acting  on  the  receptor-ligand  bonds.  C  can  be  solved  through  the  balance 
of  moment  if  the  normal  stress  distribution  on  the  adherent  membrane  is 
known  or  assumed.  Using  the  parameters  listed  in  Table  3.4,  Boltzman's 
constant  =  1 .38  x  10  ”  J/K,  and  T  =  310  K,  the  parameter  U,  representing 
the  strength  of  the  bond  responsible  for  adhesion  was  obtained  for  each  of 
the  three  treatments  with  significant  In(tp)  versus  regressions;  control, 
FMLP,  and  IL-1.  Bell  (1978)  interpreted  the  kinetic  theory  of  fracture,  such 
that,  at  r  equal  to  r,,  U,  equalled  the  product  of  the  maximum 
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force  per  bond  (fg)  and  the  maximum  stretch  per  bond  (rj.  Therefore, 

.  (3.6) 

Using  the  values  of  U„  calculated  above  and  r^  =  0.5  x  10'^  cm,  f„  for 
control,  FMLP,  IL-1,  CB,  and  colchicine  was  calculated.  Table  3.5  lists  the 
calculated  mean  values  for  each  of  the  three  treatments,  as  well  as,  the 
standard  deviations.  Figure  3.7  compares  the  values  graphically.  The  values 
of  both  U„  and  f^  were  statistically  significantly  greater  for  FMLP  than  for 
either  IL-1  or  control.  However,  the  increase  in  and  f^  was  only 


Treatment 

u. 

f. 

(dyne-cm) 

(dyne/bond) 

Control 

8.192  X  10  ± 

1.638  X  10  ®  ± 

n  =  18 

0.060  X  10  ’®SD 

0.012  X  10  ®  SD 

FMLP 

8.417  X  10  ’*  ± 

1.683  X  10  ®  ± 

n  =  13 

0.035  X  10  ’*  SD 

0.007  X  10  ®  SD 

IL-1 

8.231  X  10’*  ± 

1.646  X  10  ®  ± 

n  =  12 

0.043  X  10  ’*  SD 

0.009  X  10  ®  SD 

CB 

8.154  X  10  ’*  ± 

1.631  X  10  ®  ± 

n  =  13 

0.102  X  10  ’*  SD 

0.020  X  10  ®  SD 

Colchicine 

8.100  X  10  ’*  ± 

1.620  X  10  ®  ± 

n  =  11 

0.115  X  10  ’*  SD 

0.023  X  10  ®  SD 

aiues  Shown  are  mean  ±  standard  deviation  (SD). 


Figure  3.7  (a)  Adhesion  energy  and  (b)  bond  force  for  control. 
FMLP,  IL-1 ,  CB,  arxf  colchicine  treatment.  and  f^  were 
significantly  greater  for  FMLP  than  all  others.  Values  shown 
are  means  +/-  standard  deviations. 
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3%  greater  for  FMLP  versus  control  or  IL-1.  The  difference  in  magnitude 
was  not  significant  enough  to  conclude  whether  adhesion  through  FMLP 
stimulation  occurred  due  to  a  different  bond  or  bonding  mechanism. 
However,  the  similarity  of  Uq  and  f^  for  control  and  IL-1  stimulated  adhesion 
did  indicate  that  these  two  conditions  induced  WBC  adhesion  through  similar 
mechanisms.  U„  and  f^  were  significantly  less  for  CB  and  colchicine  as 
compared  to  the  other  three  conditions.  The  decrease  in  bond  energy  and 
bond  strength  may  by  due  to  the  WBC  and  EC  cytoskeletal  deterioration 
caused  by  the  two  treatments,  in  all  cases,  U„  and  f„  are  comparable  to 
values  calculated  by  Bell  (1977)  using  kinetic  theory. 

The  force  per  bond,  f„,  as  well  as,  the  mean  hemodynamic  values  for 
each  of  the  treatments  were  incorporated  into  a  mechanical  analysis  of  the 
force  equilibrium  required  to  keep  a  WBC  adherent  to  the  EC.  Thereby,  the 
minimum  bond  density  was  elucidated  for  control  conditions  and  FMLP  and 
IL-1  stimulated  adhesion.  Such  an  analysis  is  addressed  in  Chapter  4. 

3.2  WBC  Membrane  Mechanics 

As  a  measure  of  the  mechanical  properties  of  the  WBC  membrane,  the 
contact  length  and  angle  of  contact  between  the  WBC  and  EC  were  studied 


versus 
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Le  was  calculated  at  tp  for  each  cell  using  the  linear  regressions 
employed  earlier  to  determine  tp.  The  contact  length  at  tp  (L^tp)  was  then 
plotted  against  the  corresponding  wall  shear  stress.  Figure  3.8  and  3.9 
depict  L„  ^  versus  t^,  the  associated  first  order  linear  regressions,  and  95% 
confidence  intervals.  The  parameters,  t-ratios,  p-vaiues,  and  correlation  are 
listed  in  Table  3.6  for  each  treatment.  Although  none  of  the  regressions 
were  statistically  significant,  the  general  trends  could  be  loosely  interpreted. 
Under  control  conditions,  L^tp  increased  with  indicating  that  the  cells  were 
easily  deformable.  With  FMLP  stimulated  adhesion,  there  was  the  opposite 
trend  which  was  interpreted  as  the  ceils  that  adhered  at  higher  shear 
stresses  were  less  deformable  due  to  activation  and  subsequent 
polymerization  of  the  f-actin.  The  WBCs  under  IL*1  stimulated  adhesion 
showed  a  trend  similar  to  that  of  control.  L,,,p  increased  with  for  both  CB 
and  colchicine.  This  was  expected  since  both  of  these  treatments  have  been 
shown  to  increase  cell  deformability.  The  rate  of  increase  was  almost  four 
times  greater  for  colchicine  than  CB.  Statistical  insignificance  of  all  of  these 
regressions  was  probably  due  to  the  limited  range  of  data. 

The  contact  angle  between  WBC  and  EC  was  measured  for  all  five 
conditions.  Figure  3.10  shows  representative  trends  of  contact  angle,  B, 
versus  time  for  control,  FMLP,  and  IL-1.  Figure  3.1 1  shows  representative 
trends  for  CB  and  colchicine.  All  conditions  except  FMLP  stimulated 


Regression  parameters  and  statistics  for  ^  versus  r, 


Contact  angle  (degrees) 


Figure  3. 11  Representative  trends  of  contact  angle  versus  time 
for  WBCs  under  (a)  CB  and  (b)  colchicine  treatment. 
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adhesion  exhibited  an  exponential  decrease  in  0  versus  time  that  was  fitted 
with  the  equation: 

0{t)  =  ai^ +b  (3.7) 

where  a  and  b  are  parameters  such  that  9^,  the  initial  contact  angle,  equals 
ad  +  b)  and  a  is  the  time  constant  for  0  ^  t  ^  t,^.  The  contact  angle 
versus  time  for  FMLP  stimulated  adhesion  demonstrated  an  initial  decrease 
followed  by  a  steady  increase  in  0  to  90°  or  a  perfectly  spherical  cell.  A 
second  order  linear  regression  was  used  to  curve  fit  the  FMLP  stimulated 
WBC  change  in  0  with  time.  The  steady  increase  in  0  with  time  for  FMLP 
cells  corresponded  with  the  increased  polymerization  of  f-actin  as  the  cells 
became  more  activated.  Half  of  the  FMLP  stimulated  WBCs  exhibited  an 
initial  exponential  decrease  in  0  with  time  and  were  also  fitted  with  equation 
(3.7). 

Using  the  equations  for  0(t)  for  the  adherent  WBCs,  the  mean  value 
of  0  at  tp  was  calculated  for  each  of  the  conditions.  It  was  also  of  interest 
to  determine  if  the  cells  had  reached  a  steady  state  value  of  0  by  tp.  This 
was  determined  by  calculating  the  percent  of  total  change  in  0  that  had 
occurred  by  tp.  The  percent  of  total  change  by  tp  was  not  a  meaningful 
parameter  for  the  FMLP  cells  since  they  exhibited  an  initial  decrease  and 
subsequent  increase  in  contact  angle  with  time.  Table  3.7  lists  the  mean 
and  standard  deviation  for  the  wall  shear  stress,  contact  angle  at  tp,  as  well 
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as,  %  angle  change  by  tp,  time  constant,  a,  and  0^  for  each  condition.  The 
number  of  cells  in  each  condition  group  is  smaller  than  that  used  for  the 
contact  length  data;  this  was  due  to  the  contact  angle  not  being  discernable 
for  several  ceils  due  to  poor  image  quality.  As  seen  in  Table  3.7,  greater 
than  95%  of  total  angle  change  had  occurred  by  tp  for  WBCs  under  control 
conditions  and  IL-1  stimulated  adhesion.  Therefore,  a  steady  state  in  cell 
deformation  had  essentially  been  reached. 

Figures  3.12  and  3.13  depict  the  contact  angle  calculated  at  tp  {0^^) 
for  each  ceil  versus  the  corresponding  wall  shear  stress  for  control,  FMLP, 
IL-1,  CB,  and  colchicine.  The  parameters  and  statistical  analysis  for  each 
regression  are  listed  in  Table  3.8.  The  regression  for  control  conditions  was 
significant  at  p  <  0.02;  the  regression  for  FMLP  stimulated  adhesion  was 
significant  at  a  level  of  p  <  0.053;  the  regression  for  colchicine  treatment 
was  significant  only  at  a  level  of  p  <  0.10,  while  the  regressions  for  IL-1 
stimulated  adhesion  and  CB  treatment  were  not  significant. 

The  results  of  contact  angle  versus  wall  shear  stress  were  consistent 
with  the  earlier  conclusions  made  from  the  contact  length  data  for  control 
conditions,  FMLP  stimulated  adhesion,  and  colchicine  treatment.  WBCs 
under  control  conditions  were  readily  deformable,  WBCs  with  FMLP 
stimulated  adhesion  tended  to  adhere  at  higher  shear  rates  and  were  more 
activated,  therefore,  less  deformable,  and  WBCs  under  colchicine  treatment 
were  highly  deformable.  The  contact  angle  and  contact  length  results  for 


(b)  colchicine,  N  -  9 


Figure  3.1 3  Contact  angle  at  t  versus  wall  shear  stress  for 
(a)  CB  and  (b)  colchicine.  The  regression  for  CB  was  not 
statistically  significant,  wNIe  the  regression  for  colchicine 
was  significant  at  p  <0.11. 
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IL-1  and  CB  appeared  to  contradict  one  another.  The  lack  of  decrease  in 
contact  angle  at  tp  for  IL-1  stimulated  WBCs  with  wall  shear  stress  seemed 
to  indicate  that  the  ceils  were  less  deformable,  while  the  contact  length  data 
suggested  the  WBCs  were  readily  deformable.  The  contact  angle  data  was 
consistent  with  earlier  studies  by  Lipowsky  et  a/.  (1991).  For  CB,  the 
contact  length  data  increased  with  wail  shear  stress  which  was  congruous 
with  increased  deformability  and  previous  studies  published  by  Lipowsky  et 
al.  (1991 )  while  the  contact  angle  data  suggested  the  opposite  trend.  Since 
neither  the  regression  of  L^  ,p  versus  or  the  regression  of  0^,  versus 
were  statistically  significant  for  either  IL-1  or  CB  it  was  not  reasonable  to 
make  any  firm  conclusions. 

It  was  also  of  interest  to  investigate  how  contact  angle  calculated  at 
tp  varied  with  tp.  It  had  been  shown  in  the  present  study  that  tp  varied 
inversely  with  for  control,  FMLP  and  IL-1  and  was  virtually  Invariant  for 
colchicine  and  CB.  Contact  angle  decreased  with  increasing  for  control 
conditions  and  colchicine  treatment,  increased  with  for  FMLP  stimulated 
adhesion,  and  did  not  vary  significantly  with  for  IL-1  stimulated  adhesion 
and  CB  treatment.  Therefore,  it  was  expected  that  6  would  vary  with  tp  in 
a  direction  opposite  to  that  of  0  and  r^.  As  seen  by  Figures  3.14  and  3.15, 
the  expected  outcome  held  true  for  all  except  for  the  colchicine  treatment. 
For  control  conditions,  as  tp  increased,  indicating  a  strong  adhesive  force 
between  the  WBC  and  EC  and  increased  cell  activation,  the  contact  angle 
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also  increased.  The  regression  of  B  calculated  at  tp  versus  tp  for  control  was 
statistically  significant  for  p  <  0.05.  The  FMLP  activated  WBCs 
demonstrated  the  opposite  behavior  although  the  regression  was  not 
statistically  significant  (p  >  0.05).  The  higher  values  of  tp  were  coincident 
with  lower  shear  stresses  and  adhesion  of  less  activated  WBCs.  Since  the 
WBCs  were  not  as  activated,  they  were  more  deformabie  and  0  decreased 
as  tp  increased.  The  regression  of  9  for  IL-1  activated  adhesion  and  CB 
treatment  versus  tp  was  fairly  invariant  and  not  statistically  significant  (p  > 
0.05).  For  colchicine  treatment,  contact  angle  calculated  at  tp  decreased 
with  tp  (p  <  0.06).  This  may  indicate  that  the  cells  remained  deformable 
due  to  microtubule  disruption  regardless  of  the  degree  of  WBC  activation  and 
strength  of  adhesion  between  WBC  and  EC.  The  regression  parameters  and 
statistics  for  all  conditions  are  listed  in  Table  3.9. 

The  results  obtained  for  both  and  B  were  consistent  with  current 
knowledge  of  WBC  deformability;  thus,  both  parameters  appear  to  be  a  good 
measure  of  WBC  membrane  mechanics  and  deformability. 


Regression  parameters  and  statistics  for  versus  tp 
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CHAPTER  4 

Force  Equilibrium  Analysis  of  an  Adherent  White 
Blood  Cell  in  Shear  Flow 


Numerous  analyses  have  been  presented  in  the  literature  to  model 
bond  distribution  within  the  contact  zone  between  the  leukocyte  (WBC)  and 
endothelial  cell  membrane  during  adhesion  of  these  cells.  Most  of  these 
models  have  been  applied  to  a  leukocyte  rolling  at  a  constant  velocity. 
Dembo  et  al.  (1988)  applied  a  model  of  tape  peeling,  using  the  membrane 
equations  of  equilibrium  and  the  bond  kinetic  equations  for  adhesion  to 
determine  both  membrane  shape  and  bond  density  distribution  in  the 
adherent  region  without  hydrodynamic  effect.  Scott  (1992)  used  a  pure 
kinetic  model  of  a  rolling  WBC,  assuming  a  flat  membrane  shape  in  the 
adherent  zone.  He  estimated  a  relationship  between  the  bond  density  and 
cell  rolling  velocity.  Tdzeren  and  Ley  (1992)  used  both  bond  kinetics,  as 
well  as,  rolling  kinematics  for  their  analysis.  They  assumed  the  membrane 
shape  to  be  spherical  along  the  adherent  surface  for  low  blood  velocities. 
No  cell  membrane  equilibrium  was  involved  in  their  analysis.  The  model 
presented  here  aims  to  predict  the  adhesive  bond  distribution  which  is 
required  to  maintain  a  WBC  adherent  to  the  endothelium  based  on  the 
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mechanical  equilibrium  between  hydrodynamic  and  adhesive  forces.  This 
model,  however,  does  not  include  bond  kinetics.  The  gap  between  WBC 
membrane  and  EC  in  the  adhesion  2one  is  assumed  to  be  represented  by  an 
exponential  decay  from  the  point  of  maximum  bond  stretch  at  the  peeling 
site  to  zero  near  the  center  of  adhesion  contact. 

4. 1  Analysis 

The  analysis  presented  here  considers  the  equilibrium  of  hemodynamic 
shear  forces  versus  bond  adhesion  forces  that  must  be  maintained  in  order 
for  a  WBC  to  remain  adhered  to  the  venular  endothelium.  The  analysis  yields 
the  WBC  bond  density,  bond  normal  stress,  and  membrane  tension 
distribution  for  the  bridged  half  of  contact  length  that  is  in  tension  (Figure 
4.1).  It  is  assumed  that  the  WBC  is  stationary  (V„bc  =  0)  and  that  there  is 
no  peeling  velocity.  This  model  is  strictly  a  mechanical  analysis  and  does  not 
take  into  account  the  kinetics  of  bond  formation.  It  was  hypothesized  that 
with  appropriate  assumptions  for  the  shape  of  the  WBC  membrane  in  the 
zone  of  contact  the  results  obtained  through  the  mechanical  analysis  should 
be  close  to  that  obtained  through  kinetic  analysis  in  earlier  studies. 

At  equilibrium,  there  are  two  membrane  regions;  a  free  unbridged 
region,  and  a  bridged  region  where  the  membranes  are  held  together  with 
receptor-ligand  bonds.  The  bridged  region  or  contact  area  is  assumed  to  be 
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^contact  zone 
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Figure  4.1  Schematic  view  of  an  adhering  WBC.  (a)  Cross-sectional  view 
of  an  adhering  WBC.  AP  is  the  difference  between  the  inside  and  outside  of 
the  cell,  r,  is  the  shear  stress  acting  on  the  surface  of  the  WBC,  F,  is  the 
resultant  shear  force,  and  M  is  the  resultant  torque,  (b)  The  diagram  shows 
the  assumed  contact  area,  (modified  from  Schmid-Schdnbein,  1975) 


56 


circular  with  radius  a.  In  order  to  counter  the  torque  set  up  by  the  fluid 
which  the  bonds  are  assumed  to  be  in  tension,  beginning  with  the  trailing 
shear  force  acting  on  the  cell,  the  normal  stress  due  to  the  bonds  in  the 
adherent  region  must  be  non-uniform.  Only  the  portion  of  the  contact  region 
in  which  the  bonds  are  assumed  to  be  in  tension,  beginning  with  the  trailing 
edge  of  peeling  zone,  will  be  analyzed.  The  unbridged  region  and  the  leading 
edge  of  the  contact  zone  which  bonds  are  formed  are  not  considered  in  this 
analysis  except  for  requiring  membrane  tension  continuity  between  the 
unbridged  and  bridged  region. 

Kinematics  of  the  peeling  process  is  considered  in  the  lateral  plane 
parallel  to  the  direction  of  flow  and  the  WBC  membrane  is  considered  as  a 
two  dimensional  ring  with  the  contact  zone  described  by  the  coordinates  {s, 
6)  where  s  =  0  at  the  leading  edge  of  the  peeling  zone.  In  order  to  remain 
in  equilibrium  the  hemodynamic  and  membrane  (bond)  stresses,  as  shown 
in  Figure  4.2,  must  be  balanced.  The  membrane  tension,  T„,  acts  in  the 
plane  of  the  membrane.  A  transverse  shear,  Q,„,  acting  normal  to  the 
membrane,  results  from  bending  stresses  in  the  membrane  that  are  localized 
to  the  sharp  bending  adjacent  to  and  within  the  adherent  zone  (Evans, 
1 985).  The  attractive  normal  stress,  a„,  arises  from  the  forces  acting  on  the 
receptor-ligand  bonds  in  the  bridged  region  is  given  as  the  product  of  force 
per  bond  (f)  and  bond  density  (Nb),  i.e.  a„  =  f  •  Nb.  The  difference  in 
pressure  between  the  inside  and  outside  of  the  cell,  Ap,  is  assumed  to  be 
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surface,  the  attractive  stress,  a„,  which  represents  the  adhesive  bond  stresses,  the  pressure  difference,  Ap,  between 
the  inside  and  outside  of  the  ceil,  the  resultant  shear  stress,  r„  imparted  by  the  adhesion  to  balance  the  fluid  shear 
flow,  and  the  internal  shear  stress,  r.,  due  to  cell  deformation. 
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constant  throughout.  The  resultant  shear  stress,  r.,  imparted  to  the  WBC 
by  the  adhesion,  balances  the  drag  forces  exerted  by  the  blood  flow. 
The  internal  cytoplasmic  shear  stress,  Tp,  is  due  to  cell  deformation.  Both  r, 
and  Tp  are  assumed  to  be  constant  in  the  region  they  act. 

The  resultant  shear  stress  in  the  model,  r„  is  assumed  to  be  localized 
at  the  region  of  peeling  of  the  bridged  zone  due  to  the  fact  that  the  dominant 
factor  resisting  the  forward  roiling  motion  of  the  cell  occurs  mostly  at  this 
peeling  edge  (Dembo  et  a!.,  1988).  Physically,  the  bonds  at  the  peeling 
region  will  provide  the  greatest  stress  component  in  the  tangential  direction 
opposing  the  fluid  shear.  For  this  model,  the  bond  stress  is  assumed  to  act 
normal  to  the  endothelial  surface;  while,  r,  is  a  constant  stress  applied  over 
a  small  peeling  zone  of  the  adherent  region.  To  analyze  the  contributions  of 
the  peeling  process  to  static  equilibrium  of  the  entire  bridged  zone,  the  zone 
of  peeling  is  assumed  to  be  localized  over  a  distance,  b.  The  ratio  of  the 
zone  of  peeling  over  the  radius  of  contact  (b/a)  is  denoted  the  non- 
dimensional  parameter,  6,  and  is  arbitrarily  taken  between  0.01  and  0.75. 
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4.1.1  Peeling  Region 


First,  the  peeling  region  of  the  bridged  area  will  be  considered'.  The 
local  mechanical  equilibrium  of  the  membrane  is  given  by  the  following 
equations  (Evans  and  Skaiak,  1980): 

The  balance  of  forces  tangent  to  the  surface  is  given  by. 


ds 


(4.1) 


and  the  balance  of  forces  normal  to  the  surface  is  related  to  the  membrane 
tension  and  bending  modulus  by: 


(4.2) 


where  K„,  is  the  local  curvature  of  the  meridional  arc, 


d9 


(4.3) 


'  Superscript  1  will  be  used  to  denote  the  peeling  edge;  superscript  2 
will  denote  the  remaining  bridged  region. 
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The  local  balance  of  moments  for  the  membrane  surface  yields  the  relation 
that  the  transverse  shear  is  equal  to  the  gradient  of  the  curvature  multiplied 
by  the  bending  modulus  (Evans  and  Skalak,  1980). 

Q  (4.4) 

ds 


Since,  it  is  a  matter  of  debate  whether  or  not  the  WBC  membrane  can  offer 
a  significant  bending  moment,  zero  bending  modulus  will  be  assumed. 


Since,  the  curvature  of  the  bridged  zone  is  very  small, 

cr”’  -  0  s  0  •  (4.5) 


Therefore,  equations  (4.1)  and  (4.2)  become,  respectively. 


and 


as 


(4.7) 


One  method  of  solving  the  equations  is  to  assume  a  membrane  shape 
in  the  adherent  region.  Evans  (1985)  analyzed  a  similar  problem  without 
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shear  flow  and  was  able  to  solve  for  the  membrane  displacement  from  the 
equilibrium  position  with  the  approximation  that  membrane  tension  remained 
constant.  A  modified  version  of  Evans'  displacement  will  be  assumed  here; 

f  =  (4.8) 

where 

C  =  L  (4.9) 

ro 


and 


f  =  £  .  (4.10) 

a 

f  denotes  the  non-dimensionaiized  membrane  displacement  from  equilibrium 
or,  equivalently,  the  non-dimensionaiized  bond  stretch  from  zero,  as  shown 
in  Figure  4.3.  Bond  stretch  is  denoted  r;  the  maximum  bond  stretch  is  r„;  c 
is  a  non-dimensional  parameter  chosen  to  describe  the  exponential  rate  of 
decay  of  the  membrane  to  equilibrium  (equilibrium  is  assume  to  be  a  flat 
membrane  in  the  same  plane  as  the  EC);  and  f  is  the  spatial  coordinate,  s, 
in  the  meridional  plane  measured  from  the  trailing  edge  normalized  with 
respect  to  a,  the  radius  of  contact.  Substituting  the  non-dimensionaiized 
parameters  into  equations  (4.6)  and  (4.7)  yields: 
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(4.11) 


and 


■jiw  dO 

adi 


=  -  Lp 


(4.12) 


For  membrane  angles,  measured  relative  to  the  equilibrium  contact 
plane,  that  are  ^  30°,  the  angle  and  curvature  can  be  approximated  by  the 
first  and  second  spatial  derivative  of  the  displacement  (Evans,  1985). 


_de_  ^  dH 

ds  *  aVr 


(4.13) 


Therefore,  equation  (4.12)  becomes; 


^  f  0*^ 


(4.14) 


Equations  (4.11)  and  (4.14)  can  be  solved  In  terms  of  Tj”,  aj”,  r„ 
and  fp  by  imposing  the  boundary  conditions;  f  =  0,  Tj”  =  T„  and  = 
Op.  Applying  the  boundary  conditions  to  equation  (4. 1 4),  yields  an  equation 
for  Ap. 


Ap 


-  Oo 
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(4.15) 


Integrating  equation  (4.1 1 ),  the  balance  of  forces  tangent  to  the  membrane, 
yields: 

=  .  (4.161 


Substituting  for  Ap,  the  baiance  of  forces  normai  to  the  membrane  becomes. 


a  VJ*  -aj 


+  • 


(4.17) 


Setting  equation  (4.16)  and  (4.17)  equal,  the  rrj’*  distribution  can  be  solved 
for  the  peeling  region  of  the  bridged  zone. 


=  ((r.-r^)af  +  T^^  -e'f))  + 

a 


(4.18) 


At  this  point,  a  linearized  relation  is  introduced  between  the  molecular 
adhesion  force  and  membrane  displacement.  The  relationship  holds  for  the 
entire  bonded  region.  Such  that. 


f=fo'^  0<f=s1 


(4.19) 
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and 

f=0  r  >  1  (4.20) 

where  f  is  the  force  per  bond,  and  is  the  maximum  force  per  bond. 

Thus,  the  attractive  bond  stress,  a„,  is  equal  to  the  product  of  the  force  per 
bond  and  the  bond  density  (Nb(s)). 

a„  =  /Vb(f)  •  f  =  Mbit)  •  fo^  =  Nb{s)f„e-‘f  (4.21) 

Now,  Nb(s)'^'  can  be  solved  for  the  peeling  region. 

/V6(f)<’»=f^  ((r.-r,)af  +  (4.22)  — 


4.1.2  Bonded  Region  -  Outside  of  the  Peeling  Region 

The  assumptions  made  in  the  analysis  of  the  peeling  region  concerning 
bending  modulus  and  membrane  curvature  hold  true  for  this  region.  Hence, 
the  simplified  equations  of  static  equilibrium  for  this  region  are  as  follows: 

First,  the  membrane  equation  in  the  tangential  direction, 


adi  " 
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(4.23) 


Next,  the  membrane  equation  in  the  normal  direction, 


(4.24) 


In  order  to  ensure  continuity  between  the  peeling  region  and  the  rest 
of  the  bridged  membrane,  the  boundary  condition;  at  f  =  b/a  =  S,  Tj^(6) 
3=  is  imposed.  Since  Ap  is  assumed  to  be  constant,  it  is  the  same 

In  this  region  as  in  the  peeling  region.  Therefore,  substituting  in  for  Ap  and 
applying  the  boundary  conditions,  equations  (4.23)  and  (4.24)  become, 

r”=r^a{S-^+rji'{6)  (4.25) 


and 


'  m 


<7’*^  +  7*  ^  r  -o 

n  9  •%  o  o 

a* 


_r  e 
a*  " 


-ef 


(4.26) 
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Setting  equations  (4.25)  and  (4.26)  equal,  and  Nb'^  can  be 
solved. 


a 


a-‘ 


(4.27) 


/vb'« = I  rfl(*y  -  f)  +  rjrf)  ] 
a^fo 


(4.28) 


is  solved  using  the  momentum  balance  between  the  hemodynamic 
forces  and  the  sum  of  the  forces  due  to  the  tensile  and  compressive  bonds 
across  the  entire  contact  area.  The  stress  distributions  for  the  tensile  and 
compressive  bonds  are  assumed  to  be  equal  but  opposite.  The  moment  is 
taken  about  the  center  of  the  contact  area  so  that  the  moment  arms  are 
equal.  Therefore, 

F,' r,  *  M  =  ll  a„  X  dA  (4.29) 

A 


where  F,  is  the  force  exerted  on  the  cell  via  the  shear  flow,  r,  is  the  radius 
of  the  cell,  and  M  is  the  resultant  torque.  The  coordinate  x  is  taken  from  the 
center  of  the  contact  area.  In  order  to  solve  the  momentum  balance,  F,  and 
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M  are  approximated  by  the  following  equations  taken  from  work  by  Tozeren 
and  Skalak; 

F.  =  C,6;rr^rf  (4.30) 

and 

M  =  CiAnry^  (4.31) 

where  C,  and  Cj  are  constants.  To  simplify  the  calculation  of  equation 
(4.29),  the  approximation  f  «»  1  -  x/a  can  be  made.  By  employing  a 
trapezoidal  numerical  approximation,  equation  (4.29)  is  solved  for  f  =  0  to 
1  where  =  cr„*^'  for  0  ^  f  ^  <5  and  a„  =  for  <5  s  s  1 ,  thereby, 
elucidating  the  value  of 

4.2  Model  Validation 

The  membrane  tension,  attractive  normal  stress,  and  bond  density 
distribution  are  calculated  with  a  Fortran  program,  DIST.EXE  (Appendix  C), 
using  the  parameters  listed  in  Table  4.1 .  Nb„,  the  maximum  available  bond 
density,  is  included  In  the  model  so  that  the  bond  density  is  capped  within 
a  physiologic  range.  If  a  greater  bond  density  is  required,  the  cell  will  no 


Table  4.1  69 

Parameter  estimates 


1  Parameter 

Definition 

Range 

Reference 

ro 

maximum  bond 

interaction 

length 

5.0  X  10®  cm 

Bell,  1978 

fo 

maximum  bond 
force 

1.2  X  10® 
dyne/bond 

Bell,  1978 

•■c 

WBC  radius 

4.0  //m 

Schmid- 

Sch5nbein, 

1975 

To 

membrane  tension 
outside  bridge 
region 

0.031  -  0.035 
dyne/cm 

Dong,  1988 
Evans,  1989 

wall  shear  stress 

0.8  -  6.0 
dyne/cm* 

current 

study 

shear  stress  due 
to  cell 
deformation 

0-3.0 

dyne/cm* 

current 

study 

1  ^ 

radius  of  WBC- 
endothelium 
contact  area 

2.0  -  6.0  ^m 

current 

study 

6 

peeling  region  (b/a) 

0.01  -  0.75 

current 

study 

c 

rate  of  decay  of 
membrane  shape 

4-  7 

current 

study 

c, 

constant  that 
depends 

on  the  cell  shape 

1.7005 

Tozeren, 

1977 

C, 

constant  that 
depends 

of  the  cell  shape 

0.944 

Tozeren, 

1977 

Nb^ 

maximum  available 
bond  density 

10’° 

Dembo,  1988 

longer  be  in  equilibrium.  If  the  bond  density  distribution  does  not  change  to 
accommodate  the  shortfall,  the  cell  will  be  swept  away  by  the  shear  flow. 

Since  the  analysis  contains  two  independent  variables,  6  and  c,  it  is 
of  interest  to  ascertain  their  effect  on  the  model.  The  following  set  of 
parametric  curves  were  run  keeping  all  input  parameters  constant  while 
varying  only  the  independent  parameter  of  interest.  Figure  4.4  shows  the 
membrane  displacement  from  equilibrium  for  various  values  of  parameter  c. 
Physically,  the  parameter  c  will  be  dependent  on  both  hemodynamic  forces 
and  cell  deformability.  The  more  deformable  the  cell  the  flatter  the 
membrane  will  tend  to  be;  the  higher  the  value  of  c.  As  seen  in  Chapter  3, 
WBC  deformability  is  influe(>  ed  ‘'•j  shear  stress,  as  well  as,  chemical 
treatment  and  degree  of  cell  activation.  However,  the  stronger  the  fluid 
shear  force  the  more  the  WBC  membrane  wiil  be  peeled  away  from  the  EC 
and  the  lower  the  vaiue  of  c.  As  shown  by  the  figure,  as  c  increases,  the 
rate  of  exponential  decay  increases  and  the  final  displacement  at  f  =  1 .0  is 
closer  to  zero.  Figure  4.5  is  the  bond  density  model  sensitivity  to  the  value 
of  c.  As  c  increases,  the  bond  stretch  associated  with  membrane 
displacement  becomes  less  as  f  approaches  1 .  Since  the  force  per  bond  has 
been  linearized  with  the  degree  of  stretch,  the  force  per  bond  also  decreases; 
therefore,  the  bond  density  increases  with  c  to  compensate.  Figure  4.6 
shows  the  sensitivity  of  a„  to  parameter  c.  As  c  increases,  the  value  of  a, 
increases,  while  the  final  value  at  f  =  1  decreases.  When  c  increases,  the 
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effective  area  in  which  the  bonds  providing  the  majority  of  the  force  required 
to  counter  the  torque  set  up  by  the  shear  flow  decreases.  Since  the  required 
force  is  not  changing,  the  initial  stress  increases.  At  c  =  6,  as  shown  in  the 
figure,  u„  severely  drops  off  just  before  f  =  1 ;  this  is  due  to  the  physiologic 
bond  density  maximum  being  reached.  The  bond  density  at  this  point  can 
not  increase  any  further  to  compensate  for  the  decreasing  value  of  force  per 
bond.  The  membrane  tension  is  not  effectively  changed  by  the  value  of  c 
since  it  is  dependent  only  on  r,,  r^,  and  T^. 

Variation  in  the  independent  parameter  S  does  not  effectively  change 
the  calculation  of  Nb  (Nb  decreases  by  0.012%  from  analysis  when  S  = 
0.01  versus  6  =  0.75).  The  model  of  Is  also  fairly  insensitive  to  changes 
in  6  since  the  bond  normal  stress  is  just  the  product  of  bond  density  and 
force  per  bond.  Figure  4.7  shows  the  model  membrane  tension  sensitivity 
to  6.  The  smaller  the  area  that  r,  acts,  the  quicker  the  membrane  tension 
reaches  Its  maximum  value.  The  maximum  value  of  membrane  tension 
increases  linearly  as  6  increases.  This  trend  is  the  direct  result  of  the 
membrane  equation  governing  the  membrane  tension  in  the  tangential 
direction. 

Figure  4.8  through  4.11  show  the  membrane  displacements,  bond 
density,  attractive  normal  stress,  and  tension  distribution  for  a  series  of 
shear  stresses  (1,  2,  4,  and  6  dyne/cm*)  assuming  the  same  membrane 
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Figure  4.8  Membrane  distributions  for  ,.^  =  1.0  dyne/cm  and  zero  internal  shear  stress.  Analysis 
was  performed  with  f  =  1.2  x  10  *  dyne/bond  and  =  0.035  dyne/cm. 
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Figure  4.9  Membrane  distributions  for  s  2.0  dyne/cm  and  zero  internal  shear  stress.  Analysis 
was  performed  with  f  =  1 .2  x  10  ^  dyne/bond  and  T  s  0.035  dyne/cm. 
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Figure  4.10  Membrane  distributions  for  =  4.0  dyne/cm  and  zero  internal  shear  stress.  Analysis 
was  performed  with  f  =  1.2  x  10  dyne/bond  and  T  s  0.035  dyne/cm. 


79 


2  (uio/euAp)  “i  w 

o  o 

•  • 


rx 

<0 

to 

"tt 

CO 

0) 

0) 

0) 

0) 

0) 

<M 

CM 

CM 

CM 

CM 

(  ujo/ouAp) 

3 


CO 

■« 

s 

< 

s 

s 

c5 

o 

£ 

CO 

a 

E  . 
S  E 

C 

2  c 

o  > 

N  "O 
*5  lO 
-8 

-2  >■ 

■O  "O 

O  § 

CO  ^ 

n  o 

4 

£  -o 

CO  10 
12 
1  X 
== 

«  - 
•5  n 
c  o 
i  £ 

2*0 

« 

^  E 
»“  o 

i 

e  S. 

3  CO 

OI  « 

cE  ^ 


80 

displacement  function,  contact  radius,  and  peeling  edge,  is  zero  in  all 
fourcases.  The  membrane  displacement  was  chosen  so  the  bonds  at  the 
center  of  the  cell  contact  area  were  virtually  unstretched;  therefore, 
providing  very  little  force  per  bond.  As  discussed  earlier,  wall  shear  stress 
will  influence  the  membrane  displacement;  however,  for  this  comparison  the 
value  of  c  was  kept  constant  for  each  trial. 

As  shear  stress  increases,  so  does  the  number  of  bonds  required  to 
maintain  the  cell  in  equilibrium.  The  bond  distribution  in  each  case  rises 
exponentially  from  nearly  zero  at  the  peeling  edge  of  the  bridged  membrane 
(f  =  0)  to  a  maximum  at  the  center  of  the  contact  area.  For  all  cases,  the 
bond  density  at  the  center  of  the  contact  area  is  less  than  Nb„;  hence,  the 
estimates  given  by  the  analysis  are  within  physiologic  conditions.  The 
normal  stress  distributions  have  the  same  profiles  but  different  magnitudes 
for  all  shear  stresses;  a„,  the  stress  at  f  =  0  increases  proportionally  with 
the  level  of  wall  shear  stress.  The  normal  stress  distribution  is  similar  to 
thatpredicted  by  Hammer  and  Lauffenburger  (1989).  Membrane  tension 
increases  slightly  within  the  peeling  zone  and  then  remains  constant  for  the 
remainder  of  the  bridged  region. 

The  conditions  in  Figure  4.1 1  were  duplicated  with  the  addition  of  a 
non-zero  in  the  direction  opposite  of  r,  to  take  Into  account  a  internal 
cytoplasmic  shear  stress  caused  by  cell  deformation  (Figure  4.12).  The 
attractive  stress  distribution  and  bond  density  distributions  do  not  change 
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Figure  4.12  Membrane  distributions  for  r..  =  6.0  dyne/cm  and  r..  =■  3.0  dyne/cm  .  Analysis 

-5  P 

was  performed  with  f  =  1 .2  x  10  dyne/bond  and  T  =  0.035  dyne/cm. 
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significantly;  however,  there  is  a  steady  decay  in  membrane  tension  as  f 
approaches  1  outside  the  immediate  peeling  zone  ^  f  ^  1).  These 
results  appear  to  indicate  that  as  the  cell  deforms  and  an  internal  shear 
stress  on  the  cell  membrane  within  the  contact  region  is  produced  in  the 
same  direction  of  the  fluid  flow,  the  membrane  tension  decreases  from  the 
maximum  at  the  peeling  edge  to  a  minimum  at  the  leading  edge. 

In  order  to  assess  the  effect  of  deformation  on  the  bond  density 
distribution,  the  mechanical  equilibrium  equations  were  calculated  over  a 
range  of  contact  lengths.  Shear  stress  (r^)  was  held  constant  and  the 
peeling  region  was  equal  to  10%  of  the  contact  radius  for  all  cases.  Once 
again,  although  deformation  will  influence  the  membrane  displacement,  the 
parameter  c  defining  the  rate  of  exponential  decay  for  the  membrane 
displacement  was  held  constant.  The  results  are  shown  in  Figure  4.13(a}, 
the  bond  density  versus  f  and  4.13(b),  a  blow-up  of  (a)  at  the  peeling  region 
(0  f  :S  J).  As  the  contact  length  increases,  the  bond  density  required  to 
maintain  the  cell  in  equilibrium  decreases  at  the  center  of  the  contact  area, 
as  does  the  bond  density  at  the  peeling  edge.  This  result  is  consistent  with 
the  assumption  that  as  contact  area  increases  for  the  more  deformable  cells, 
the  adherent  WBC  will  be  able  to  sustain  higher  shear  forces.  Thus,  contact 
area  is  a  chief  determinant  of  adhesive  force. 

As  seen  in  Chapter  3,  the  value  of  f^,  the  force  per  bond,  can  vary 
with  chemical  treatment.  Therefore,  it  is  of  interest  to  examine  the  change 


Figure  4.13  Effect  of  cor^tect  length  on  bond  density  distribution, 
(a)  is  full  scale  and  (bl  is  a  blow-up  of  the  peeling  edge.  Analysis 
was  performed  using  the  following  parameters;  c  =  5,tf=0.10. 

^  =  6.0  dyne/cm  ,  t  =  0,  f^  =  1.2  x  lO^dyne/bond,  an 
t],  =  0.035  dyne/cm.’^ 
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in  bond  density  distribution  as  a  result  of  varying  f„  while  keeping  all  other 
parameters  constant  (Figure  4.14).  As  expected,  both  the  initial  and  final 
value  of  bond  density  decrease  as  the  force  per  bond  increases.  However, 
within  the  small  range  of  values  (1.5  ^  fo  ^  2.0  x  10'^  dyne/bond)  seen  in 
Chapter  3  the  variation  in  bond  density  is  not  great. 

It  is  now  of  interest  to  compare  the  results  obtained  through  this 
mechanical  analysis  to  results  obtained  incorporating  the  kinetic  behavior  of 
the  receptor-ligand  bonds.  As  mentioned  before,  Tdzeren  and  Ley  (1992) 
analyzed  the  bond  density  distribution  with  a  model  taking  into  consideration 
bond  kinetics  and  rolling  kinematics.  Figure  4.15  is  a  schematic  of  a  rolling 
leukocyte  used  for  their  analysis.  The  model  assumed  the  membrane  shape 
to  be  a  rigid  cylinder  and  that  the  bonds  had  the  mobility  necessary  for  bond 
formation  but  no  diffusivity.  As  the  kinematics  of  rolling  dictate,  the 
adhesion  bonds  were  first  compressed  than  stretched  during  roliirig.  The 
compressed  bonds  were  governed  by  both  an  attachment  rate  and 
detachment  rate,  while  the  stretched  bonds  were  only  governed  by  a 
detachment  rate.  Figure  4.16  shows  the  bond  density  distribution  results 
Tdzeren  and  Ley  obtained.  Implicit  in  their  model  was  the  assumption  that 
the  available  leukocyte  receptors  interacted  with  their  counter-receptors  at 
any  one  instant.  Therefore,  the  bond  density  reached  the  same  maximum 
for  both  cases.  The  maximum  bond  density  equaled  the  surface  density  of 


Figure  4.15  The  schematic  diagram  of  a  leukocyte  rolling  Figure  4. 16  The  surface  density  of  adhesion  bonds 
on  vascular  endothelium.  (Top)  Geometric  parameters  as  a  function  of  distance  along  the  vascular 

used  in  the  analysis  are  shown.  (Bottom)  Configuration  endothelium  (x).  (Reproduced  from  Tozeren  and 

of  adhesion  bonds  in  the  contact  area  is  shown  schem-  Ley,  1992) 

atically.  (Reproduced  from  Tozeren  and  Ley,  1992) 
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endothelial  receptors;  the  lower  of  the  two  receptor  densities,  leukocyte 
andendothelial.  The  number  of  stretched  bonds  in  the  detachment  zone  (- 
1.2  ^  X  ^  0)  increased  with  increasing  shear  rate. 

Figure  4. 1 7  and  4. 1 8  show  the  results  obtained  by  the  mechanical 
analysis  presented  here  for  similar  conditions  as  used  by  Tdzeren  and  Ley. 
The  results  obtained  by  Tdzeren  and  Ley  are  also  shown  on  the  figure  for 
comparison.  The  value  of  c  used  in  Figure  4. 1 7  is  such  that  the  bond  length 
at  the  center  of  the  contact  area  is  very  close  to  zero.  Therefore,  the  bond 
density  rose  exponentially  from  a  non-zero  value  at  the  peeling  edge  to  the 
physiologic  maximum,  Nb„,  prior  to  the  center  of  the  contact  area.  In  order 
to  maintain  the  WBC  adherent  to  the  endothelium,  one  would  expect  that  the 
bond  distribution  would  shift  to  the  left  or  the  membrane  shape  would 
change  to  compensate  for  the  bond  density  shortfall.  Similar  to  the  results 
obtained  by  Tdzeren  and  Ley,  the  number  of  bonds  within  the  contact  region 
increases  with  shear  stress.  To  keep  the  modei  results  within  physiologic 
range,  the  bonds  at  the  center  of  the  contact  area  were  given  a  small  initial 
amount  of  stretch  by  decreasing  the  parameter  c.  It  is  obvious  from  Figure 
4.18  that  the  bond  density  required  to  keep  the  ceil  in  equiiibrium  is 
drastically  reduced  when  the  bonds  at  the  center  can  produce  a  greater 
force.  Figure  4.18  also  shows  the  same  trends  with  shear  stress  as  Figure 
4.17.  The  mechanical  model  presented  here  correlates  well  with  the  kinetic 
and  kinematic  model  introduced  by  Tdzeren  and  Ley. 
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distribution  along  the  adherent  megnbrane  for  the  parameter  values;  c  =  5,  a  =  1 .2 
a  =  0.10,  Tp  *  0,  =  1.2  X  10  dyne/bond.  and  =  0.035  dyne/cm.  Wall  shear 

stresses  and  contact  area  radius  are  the  same  as  used  by  Tozeren  and  Ley  (1992). 
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Figure  4. 1 8  Bond  density  distribution  for  larger  membrane  gap.  Bond  density 
distribution  along  the  adherent  membrane  for  the  parameter  values;  c  »  4,  a  =  1 .2  ^ 
6  =  0.10,  Tp  *  fo  =  1.2  X  10  dyne/bond,  and  =  0.035  dyne/cm.  Wall  shear 

stresses  and  contact  area  radius  are  the  same  as  used  by  Tozeren  and  Ley  (1992). 
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Scott's  pure  kinetic  model  also  produced  a  similar  bond  density  curve 
as  obtained  here  with  the  minimum  bond  concentration  located  at  the  peeling 
edge  where  r  =  r„.  The  model  assumed  by  Scott  was  of  a  flat  linear 
membrane.  The  model  In  Scott's  work  produced  a  "fuller"  curve  where  the 
exponential  rise  in  bond  density  occurred  prior  to  the  half  way  point  of  the 
contact  radius.  As  in  both  the  mechanical  model  presented  here  and  the 
model  presented  by  Tdzeren  and  Ley,  the  total  number  of  bonds  increased 
with  shear  rate. 

Dembo  et  a/.'s  model  incorporated  both  the  peeling  mechanics  and 
bond  kinetics.  Their  analysis  also  resulted  in  an  exponential  rise  in  bond 
density  from  the  peeling  edge  reaching  a  constant  value  equal  to  the  number 
of  receptors  available  for  bonding. 

The  mechanical  analysis  presented  here  provides  remarkably  similar 
results  to  that  obtained  by  several  different  kinetic  and  kinematic  models  of 
bond  density  distribution.  It  has  been  shown  the  shape  of  the  membrane  in 
the  adherent  zone,  the  contact  length,  and  the  fluid  shear  stress  all  have  a 
great  effect  on  the  number  of  bonds  required  within  the  adherent  region  to 
maintain  the  WBC  in  equilibrium.  It  is  now  of  interest  to  incorporate  the 
experimental  results  from  Chapter  3  with  the  model  presented  here  in  order 
to  analyze  to  the  force  equilibrium  required  to  keep  a  WBC  adherent  to  the 


EC  under  several  conditions. 
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Chapter  5 

Discussion  of  Bond  Density  Distributions  for  Control  Conditions 
and  FMLP  und  IL-1  Stimulated  Adhesion 

The  in  vivo  measurements  of  transients  in  contact  zone  length  yielded 
the  force  per  bond  (f,,)  for  control,  FMLP,  and  IL-1  stimulated  adhesion. 
These  parameters,  as  listed  in  Table  5.1,  were  incorporated  into  the 
mechanical  peeling  analysis  so  that  membrane  bond  density  could  be 
calculated  for  control  conditions  and  stimulated  adhesion.  With  the  addition 
of  fo,  a  kinetic  parameter  was  included  in  the  model  that  was  directly  related 
to  the  life  time  of  the  bonds.  The  internal  shear  stress,  was  assumed  to 
be  zero  for  all  analysis. 

Figures  5.1  through  5.3  depict  the  membrane  distributions  for 
displacement,  bond  density,  normal  stress,  and  membrane  tension  for 
control,  FMLP,  and  IL-1,  respectively.  The  membrane  distributions  are  over 
the  half  of  the  contact  area  that  is  in  tension  normalized  by  the  length  of  the 
contact  radius. 

The  bond  density  was  greatest  for  FMLP  stimulated  adhesion  and  least 
for  IL-1  stimulated  adhesion.  Since  the  mean  wall  shear  stress  was  greatest 
for  FMLP,  the  bond  density  required  to  maintain  the  WBC  in  force  equilibrium 
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Figure  5.1  Membrane  distributions  for  WBC  adhesion  under  control  conditions.  Analysis  was  performed 
with  =  7.13  dyne/cm^,  a  =  4.2  yum,  a  -  0.10,  c  =  5,  f^  =  1.64  x  10  dyne/bond  and 
T  3  0.035  dyne/cm.  Average  bond  density  for  control  conditions  was  3.65  bonds/^xm  . 


(a)  MEMBRANE  J  I  BOND  DENSITY 
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Figure  5.2  Membrane  distributions  for  WBC  FMLP  stimulated  adhesion.  Analysis  was  performed 
with  =  1 1.65  dyne/cm^  a  =  4.4  ^m,  s  =  0.10,  c  =  5,  f^,  =  1.68  xIO'  dynes/bond  and  ^ 
T  =  0.035  dyne/cm.  Average  bond  density  for  FMLP  stimulated  adhesion  was  5.05  bonds/^^m 
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Figure  5.3  Membrane  distributions  for  WBC  IL-1  stimulated  adhesion.  Analysis  was  performed 
with  *  7.43  dyne/cm^,  a  =  5.1  /jn.  ^  =  0.10,  c  =  5,  =  1.65  xlO  dynes/bond  and  ^ 

T  3  0.035  dyne/cm.  Average  bond  density  for  IL-1  stimulated  adhesion  was  2.12  bonds/^m  . 
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was  greatest  even  though  the  force  per  bond  was  also  the  largest  for  FMLP. 


In  order  to  accurately  assess  the  contribution  that  the  increased  value  of  f, 
played  in  determining  the  average  bond  density  for  FMLP  stimulated 
adhesion.  Figure  5.4  depicts  the  average  bond  density  versus  the  force  per 
bond,  fo,  holding  all  other  parameters  constant.  It  is  clear  from  this  figure 
that  for  f„  ^  2.5  x  10  ®  dyne/bond  small  changes  in  bond  force  will  have  a 
large  effect  on  the  average  bond  density  required  to  maintain  the  WBC  in 
force  equilibrium.  However,  for  the  very  small  variations  in  f^  determined  for 
the  three  conditions  from  the  peeling  measurements,  there  will  be  a  minimal 
effect  on  the  potential  number  of  adhesion  bonds.  IL-1  stimulated  adhesion 
and  adhesion  under  control  conditions  occurred  at  similar  hemodynamic 
forces  and  yielded  nearly  identical  values  of  force  per  bond.  Yet,  the 
required  bond  density  was  less  for  IL-1  stimulated  adhesion.  The  reason  for 
the  difference  is  due  to  the  increased  mean  contact  length  seen  with  IL-1. 
Figure  5.5  is  a  composite  of  the  bond  density  distributions  for  all  three 
conditions.  The  average  bond  densities  calculated  for  each  condition  for  a 
unit  depth  are  shown  in  Figure  5.6.  The  average  bond  density  for  FMLP  was 
5.05  bonds/|/m*  versus  3.65  bonds/^/m*  for  control  and  2.12  bonds///m*  for 
IL-1  stimulated  adhesion.  The  standard  deviations  for  the  average  bond 
densities  in  Table  5.1  and  Figure  5.6  were  derived  from  the  standard 
deviations  in  f^. 


(  uirZ/spuoq)  eOejaAV 


experimental  conditions  for  f  determined  through  the  peeling  process. 


Average  bond  density  for  control  cornlitions  and  FMLP  and  IL-1  stimulated  adhesion, 
error  bars  represent  the  variance  in  Nb  due  to  the  standard  deviation  for  f  . 
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Since  a;,  =  f  ■  Nb,  it  is  not  surprising  to  see  that  bond  normal  stress 
foiiowed  the  same  trends  as  bond  density  being  greatest  for  FMLP  and  ieast 
for  IL-I .  The  order  of  magnitude  of  normai  stress  was  within  the  same  range 
as  cited  by  Schmid-Schonbein  and  others  (1975)  for  the  normal  stress  of 
interaction  between  the  WBC  and  endotheiium.  The  totai  force  can  be 
calculated  for  each  condition  by  integrating  the  stress  distribution  over  the 
area.  Since  the  stress  distribution  was  nearly  constant,  an  approximation  of 
totai  force  was  made  by  muitipiying  the  initial  stress,  by  the  contact  area 
(ira^)  taken  at  tp.  This  yieided  the  total  force  for  each  condition  required  to 
maintain  the  WBC  in  static  equilibrium  with  the  EC  as  follows;  for  control 
conditions  -  1.13  x  10^  dynes,  for  FMLP  -  1.76  x  10’^  dynes  and  for  IL-1 
stimulated  adhesion  -  1 .539  x  10  *  dynes.  The  order  of  magnitude  obtained 
for  force  were  consistent  with  in  vitro  experimentai  data  obtained  by  Tdzeren 
et  aL  (1989)  for  cell  to  cell  adhesion  of  a  cytotoxic  T-cell  and  a  target  cell 
separated  using  micromanipulation  techniques. 

While  the  normal  stress  was  almost  twice  as  great  for  control  versus 
IL-1,  the  membrane  tension  distributions  were  comparable.  In  the  present 
analysis,  the  membrane  tension  was  dictated  solely  by  the  hemodynamic 
stress,  T„,  and  the  distance  along  the  peeling  edge,  6.  Therefore,  bond 
density  and  normal  stress  distribution  did  not  appreciably  affect  T„.  Since 
the  wall  shear  stress  for  control  and  IL-1  stimulated  adhesion  was  virtually 
the  same,  so  were  the  T„  distributions. 


A  number  of  inferences  can  be  made  from  the  discussion  above.  First, 
the  increased  strength  of  adhesion  seen  with  FMLP  does  not  appear  to  be 
due  to  the  avidity  of  the  bond  but  rather  the  increased  bond  density. 
Secondly,  deformability  plays  a  significant  role  in  decreasing  the  required 
bond  density  by  increasing  the  contact  area  between  the  WBC  and  EC  as 
clearly  shown  in  the  comparison  between  the  control  condition  and  IL-1 
stimulated  adhesion. 


Chapter  6 

Summary  and  Conclusion 
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Both  WBC'EC  contact  angle  and  contact  length  versus  wall  shear 
stress  gave  a  relatively  good  measure  of  cell  deformability  that  was 
consistent  with  previously  reported  data.  Contact  length  increased  while 
adhesive  bonds  were  being  formed  and  decreased  when  the  rate  of  bonds 
breaking  exceeded  the  rate  of  bond  formation.  This  transition  time  was 
identified  as  the  time  to  peeling  coincident  with  the  maximum  contact 
length.  The  disruption  of  WBC-EC  bonds  as  indicated  by  the  time  to  peeling 
appeared  to  be  consistent  with  the  kinetic  theory  of  fracture.  Using  the 
kinetic  theory  of  fracture,  the  bond  adhesion  energy  and  bond  force  under 
control  conditions,  FMLP  and  IL-1  stimulated  adhesion,  and  CB  and 
colchicine  treatment  were  elucidated.  The  bond  adhesion  energy  and  bond 
force  were  significantly  greater  for  FMLP  as  compared  to  control  and  IL-1; 
however,  not  by  a  great  enough  magnitude  to  conclude  whether  a  different 
type  of  bond  or  bonding  mechanism  was  involved.  Since  the  bond  adhesion 
energy  was  statistically  the  same  for  both  control  and  IL-1  activated 
adhesion,  it  was  concluded  that  both  conditions  induced  WBC  sticking 
through  similar  mechanisms. 
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A  theoretical  analysis  was  presented  in  which  the  bond  density, 
normal  stress,  and  membrane  tension  distributions  were  elucidated  for  an 
adherent  WBC  in  shear  flow.  The  WBC  membrane  was  considered  as  a 
cortical  shell  with  an  assumed  shape.  In  order  to  maintain  the  ceil  in  static 
equilibrium,  the  adhesive  forces  due  to  the  bonds  keeping  the  WBC  stuck  to 
the  EC  were  required  to  equal  the  hemodynamic  forces  tending  to  peel  the 
WBC  off  the  EC.  It  was  seen  that  both  shear  stress  magnitude  and 
deformability,  as  indicated  by  contact  length,  were  important  parameters  in 
determining  the  bond  density  distribution.  The  model  provided  an  estimate 
of  bond  density,  normal  stress,  and  membrane  tension  that  was  consistent 
with  other  kinetic  based  analysis. 

The  values  of  f^  determined  experimentally  for  control,  FMLP,  and  IL-1 
were  incorporated  into  the  model  along  with  the  mean  hemodynamic 
parameters  for  each  of  these  conditions.  The  inclusion  of  f^  into  the 
mechanical  model  incorporated  a  kinetic  parameter  that  was  directly  related 
to  the  bond  life  time.  The  average  bond  density  varied  for  control  versus  IL- 
1  even  though  the  hemodynamic  force  and  bond  force  were  nearly  identical. 
The  average  bond  density  was  lower  for  IL-1  due  to  the  increased  contact 
area.  Finally,  since  the  avidity  of  the  receptor-ligand  bonds  for  integrin 
(FMLP)  versus  selectin  (IL-1)  adhesion  was  similar  in  magnitude,  the  FMLP 
adhesion  appeared  to  be  much  stronger  due  to  the  increase  in  average  bond 
density  as  calculated  using  the  theoretical  model  rather  than  the  bond  force. 
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One  may  speculate  as  to  the  plausible  physiologic  consequence  of  this 
conclusion.  Selectin  mediated  adhesion,  as  demonstrated  with  IL-1 
treatment,  occurs  in  response  to  low  levels  of  antigens  and  primes  the 
system  for  a  greater  immune  response  through  integrin  mediated  WBC 
adhesion.  This  pre-activation  includes  the  influx  of  WBCs  into  the  area  of 
inflammation  and  increased  WCB  roiling  along  the  post-capillary  venuiar 
endothelium.  This  is  in  accord  with  the  lower  levels  of  adhesive  force  seen 
in  the  present  study  for  IL-1 .  Integrin  mediated  adhesion,  as  demonstrated 
with  FMLP,  is  the  primary  cellular  immune  response  resulting  in  not  only 
increased  WBC  adhesion  but  greater  cell  motility;  therefore,  precipitating  the 
emigration  of  WBCs  through  the  microvessel  walls.  In  this  light,  the  results 
obtained  for  FMLP  stimulated  adhesion  are  consistent  with  the  critical  role 
integrin  mediated  adhesion  plays  in  the  body's  response  to  acute  infection. 

The  present  study  logically  leads  to  ideas  for  future  research.  It  would 
be  advantageous  for  both  the  study  of  the  peeling  process  and  WBC 
deformability  to  control  the  wail  shear  stress  using  a  technique,  such  as, 
micro-occlusion.  The  use  of  monoclonal  antibodies  to  block  one  or  more  of 
the  receptors  on  the  WBC  or  EC  could  lead  to  a  more  precise  interpretation 
of  the  role  of  the  individual  adhesion  molecules.  Finally,  the  theoretical 
analysis  could  be  improved  so  that  the  membrane  shape  could  be  determined 
analytically  and  cell  deformability  included. 
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Appendix  A 
Hemodynamic  Data 


A.1  Hemodynamic  Data  •  Control 


CELL  ID  DATE  TIME  DIAM  SHEAR  RATE 
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A.2  Hemodynamic  Data  •  PMLP 
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Appendix  B 
Experimental  Analysis 


B.3  Experimental  Analysis  -  IL-1 
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Values  could  not  be  determined  due  to  image  quality. 


B.5  Experimental  Analysis  -  Colchicine 
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Appendix  C 

Bond  Distribution  Program 
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Bond  Distribution  Program 


PROGRAM  DIST 

C  PROGRAM  SOLVES  NB,  STRESS,  AND  TENSION  DISTRIBUTION  FOR 
C  WBC  EC  ADHESION. 

C  WRITTEN  BY  ERIKA  J.  STRUBLE,  29  JANUARY  1993 
C  MODIFIED  1 5  FEBRUARY  1 993  TO  INCLUDE  TWO  REGIONS 
C  . 

REAL  C.A,TW.TP,RO,FO.C1,C2,RC,PI.FS.M.TS,TO.DLTAX,X,DLTAP,F1JX. 

•  FI  JXX.F1  JXN.IF1  ,F2  JX,F2JXX.F2JXN.IF2,SIG0.S,NB,SIGN.TM,S1  ,b. 

•  FI  XB,F2XB.F3JX,F3JXX.F3JXN.IF3,TMb.SS.NBMAX.NB1  b.NBI  0,NB2A, 

•  NB2b.NNB.CN.BN,ZETA,ZETAN.AU 

INTEGER  N,J.Z,I 
CHARACTER  *20  FLNAM 
DATA  IOUT/21/ 

C  FORMAT  STATEMENTS 

7  FORMAT!'  ',F4.2,'  •,F20.4,'  'F9.4.'  '.F9.4,'  ',F9.4,'  'F9.4,'  ' 

•  .F9.4) 

C  PROGRAM  ASKS  FOR  THE  FOLLOWING  INFORMATION: 

C  FLNAM  =  FILE  TO  CONTAIN  NB.  SIGMA,  AND  TENSION  DISTRIBUTION 
C  CN  -  NON-DIMENSIONAL  PARAMETER  DESCRIBING  EXPONENTIAL  DECAY 
C  A  =  RADIUS  OF  CONTACT  AREA  (CM) 

C  AU  RADIUS  OF  CONTACT  AREA  (uM) 

C  TW  =  WALL  SHEAR  STRESS  (DYNE/CM2) 

C  TP  •=  INTERIOR  SHEAR  STRESS  (DYNE/CM2) 

C  bN  =  NON-DIMENSIONALIZED  REGION  OF  PEELING 

WRITE(*,*)'INPUT  VALUE  OF  CN:' 

READC,*)CN 

WRITEC,*)'INPUT  RADIUS  OF  CONTACT  AREA  (urn):' 

READC,*)AU 

WRITEC,*)'INPUT  VALUE  OF  bN:' 

READC,*)bN 

WRITE|*,*)'INPUT  VALUE  OF  TW:' 

READ(*,*)TW 

WRITE(*,*)'INPUT  VALUE  OF  TP:' 

READC,*)TP 

WRITEC,*)'INPUT  VALUE  OF  FO  (DYNE/BOND):' 

READ(*,*)FO 
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C  DIMENSIONALIZE  CN  AND  bN 

A  =  AU*lE-4 
C  =  CN/A 
b  =  bN*A 

C  DEFINE  CONSTANTS 

aNAM  =  'A;DIST.OUT' 
RO=0.5E-7 
Cl  =1.7005 
C2  =0.944 
RC=4.0E-4 
Pl  =  3.1415927 
TO  =0.035 
NBMAX=100 

C  CALCULATE  FS.  M.  TS 


FS  =  C1*6*PI*TW*RC**2 
M  =  C2*4*PI*TWRC**3 
TS  =  FS/(PrA**2| 

WRITEr,*)'FS  =  ',FS 
WRITEC.TM  =  ',M 
WRITEC/)'TS  =  ',TS 
WRITE|*,*)'TO  =  ',TO 

WRITEC,*)'CALCULATING  SIGO,  NB,  SIGN,  AND  TM  DISTRIBUTION' 

WRITE! •.•I'WRITING  TO  FILE  ArDIST.OUT' 

C  CALCULATE  SIGMA  0  USING  TRAPAZOIDAL  NUMERICAL  APPROXIMATION 

N  =  100 

DLTAX  =  A/100 
C  F1(X)=0 

F1XB  =  ((TS-TP)*(A-b)+T0-T0*EXP(C*(A-b)))*EXP|-C*(A-b))*b* 

*  SORT(A**2-b**2| 

C  SUMMATION  OF  FI  (J*DELTAX) 

F1JX=0 
F1JXX  =  0 


DO  10  J=t.  N-1.  1 


C  INITIALIZE  F1JXN 


F1JXN  =  0 
X=J*DLTAX 
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IF  (X  .LE.  b)  THEN 

F 1JXN  =  {(TS-TP)  *  ( A-X)  +  TO-TO  *  EXP(C  *  ( A-X)))  *  EXP(-C  *  ( A-X))  •  X« 

*  SQRT|A**2-X**2J 

F1JXX  =  F1JX  +  F1JXN 
F1JX  =  F1JXX 

ENDIF 

10  CONTINUE 

C  NUMERICAL  SOLUTION  FOR  INTEGRAL  FI  (X) 

IF1  =(DLTAX/2r(F1XB  +  2‘F1JX» 

C  F2(X) 

F2XB  =  ITS  •  b-TP*  ( A-b)  +  TO-TO  *  EXPIC*  (A-b)»  *  EXPI-C  *  I  A-b))  ’b* 

•  SORT(A**2-b**2) 

C  F2X0=0 

C  SUMMATION  OF  F2|J*DLTAX) 

F2JX  =  0 
F2JXX  =  0 


DO  20  J  =  1.N-1,1 

C  INITIALIZE  F2JXN 

F2JXN  =  0 
X=J’DLTAX 

IF  (X  .GT.  b)  THEN 


F2JXN  =  (TS*b-TP’|A-X)  +TO-TO’EXP|C’|A-X)))*EXP|-C*|A-X))’X* 
SQRT(A**2-X**2) 

F2JXX«F2JX  +  F2JXN 
F2JX  =  F2JXX 

ENDIF 


20  CONTINUE 
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C  NUMERICAL  SOLUTION  FOR  INTEGRAL  F2(X) 

IF2  =  |DLTAX/2)*|F2XB  +  2*F2JX) 

C  F3IX) 

C  F3(0|=0 
C  F3(A)  =  0 

C  SUMMATION  OF  F3(J*DLTAX) 

F3JX=0 

F3JXX=0 

DO  30  J  =  1,N-1.1 

C  INITIALIZE  F3JXN 

F3JXN  =  0 
X=J*DLTAX 

F3  JXN  =  X*  SORT!  A  *  *  2-X  *  *  2» 

F3JXX  =  F3JX  +  F3JXN 
F3JX-F3JXX 

30  CONTINUE 

C  NUMERICAL  SOLUTION  FOR  INTEGRAL  F3|X) 

IF3  =  DLTAX*F3JX 

C  SOLVE  FOR  SIGMA  0  AND  DELTA  P 

SIGO  =  HIFS*RC  +  M)/4)-C*  •2*RO*(IF1  +  IF2))/|IF3) 

WRITE!*,  TSIGM  AO  =  ',SIG0 

DLTAP=  TO*C**2*RO-SIGO 
WRITE(*,*)'DELTA  P  =  '.DLTAP 

C  CALCULATE  NB.  SIGN,  TM,  NORMALIZED  ZETA  AND  WRITE  TO  FILE 

TMb  =  <TS-TP)*b  +  TO 

OPEN(UNIT  =  lOUT.FILE  =  FLNAMJ 
C  WRITE|IOUT,6) 

Z  =  100 
S  =  0 


WRITE(*.*)Z 
DO40l  =  0.Z,1 
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IF  (S  .LE.  b)  THEN 

C  NB  IS  CALCULATED  FOR  BONOS/UM2 

NB  =  «C**2*RO/FO)*((TS-TP)*S  +  TO-TO*EXP(C*S))  +  SIGO/(FO‘EXP 

*  (-C*S)))*1E-8 

IF  (NB  .GE.  NBMAX)  THEN 
NB=NBMAX 
ENDIF 

SIGN  =  |NB*1  E8)*FO*EXP(-C*S) 

TM  =  (TS-TP)*S  +  TO 

ZETA  =  RO*EXP(-C’S) 

ZETAN  =  ZETA/RO 

ELSE 

NB  =  l|C**2*RO/FO)*(TP*|b-S»  +  TMbHTO*C**2*RO-SIGO)/ 

•  (FO*EXP|-C*S))»*t’ ? 

IF  (NB  .GE.  NBMAXI  THEN 
NB  =  NBMAX 
ENDIF 

SIGN  =  (NB*  1 E8)  *FO*EXP(-C*S) 

TM=TP*(b-S»+TMb 

ZETA  =  RO*EXP(-C*S) 

ZETAN  =  ZETA/RO 

ENDIF 

SS  =  S*1E4 

WRITE(IOUT.7)SS.NB,SIGN,TM,DLTAP,SIGO,ZETAN 

S1=S  +  6E-6 
S-S1 
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40  CONTINUE 

CLOSE(UNIT  =  IOUT) 

STOP 

END 


